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We give updated constraints on hypothetical light bosons with a two-photon coupling such as 
axions or axion-like particles (ALPs). We focus on masses and lifetimes where decays happen near 
big bang nucleosynthesis (BBN), thus altering the baryon-to-photon ratio and number of relativistic 
degrees of freedom between the BBN epoch and the cosmic microwave background (CMB) last 
scattering epoch, in particular such that N™ B < N BBti and r/ CMB < ?; BBN . New constraints 
presented here come from Planck measurements of the CMB power spectrum combined with the 
latest inferences of primordial 4 He and D/H abundances. We find that a previously allowed region 
in parameter space near m = 1 MeV and r = 100 ms, consistent with a QCD axion arising from 
a symmetry breaking near the electroweak scale, is now ruled out at > 3<r by the combination of 
CMB+D/H measurements if only ALPs and three thermalized neutrino species contribute to N e s- 
The bound relaxes if there are additional light degrees of freedom present which, in this scenario, 
have their contribution limited to A N e s = 1.1 ± 0.3. We give forecasts showing that a number of 
experiments are expected to reach the sensitivity needed to further test this region, such as Stage-IV 
CMB and SUPER-KEKB, the latter a direct test insensitive to any extra degrees of freedom. 


I. INTRODUCTION 

Recent improvements in both measurements of the 
Cosmic Microwave Background (CMB) angular power 
spectrum and inferences of primordial elemental abun¬ 
dances formed during Big Bang Nucleosynthesis (BBN) 
motivate us to reconsider bounds on hypothetical sce¬ 
narios which alter the expansion rate or inject energy 
into the plasma around these two epochs. One of the 
simplest scenarios involves a radiatively decaying parti¬ 
cle, for which axions or so called “axion-like particles” 
(ALPs) provide a theoretically well motivated candidate. 
Axions arise from perhaps the most elegant solution to 
the strong CP problem as the pseudo Nambu-Goldstone 
boson of a new spontaneously broken symmetry [1-4]. 
Axions have a mass and standard model couplings con¬ 
trolled by a single parameter: the energy scale of the 
symmetry breaking. ALPs form a more general class of 
particles where the mass and couplings are independent. 
Such models can arise from other new symmetries which 
are spontaneously broken [5, 6], and in string theory [7]. 

In the parameter space of interest here, ALPs are 
weakly interacting, making them difficult to detect in the 
laboratory. Cosmological constraints serve as a natural 
complement as the weak coupling generally leads to later 
decays allowing the particles to become non-relativistic 
and pick up energy compared to the plasma, leading to 
observable consequences upon their decay. Astrophysi- 
cal bounds are also important as ALPs provide a new 
method for energy release from stars and supernovae. 
Some early calculations and compilation of cosmologi¬ 
cal, astrophysical, and laboratory bounds on ALPs in¬ 
clude those from Masso and Toldra [8, 9]. Cosmological 
bounds were recently updated by Cadamuro ef al. [10] 
which considered only axions, and Cadamuro and Re¬ 
dondo [11] who extended this more generally to ALPs. 


Among other significant advances, Cadamuro and Re¬ 
dondo [11] used newer data, treated out-of-equilibrium 
decays more carefully, and performed precise calculations 
of the implications for BBN. These and other known 
ALP bounds are tabulated in [12-14]. Our work updates 
these by 1) using the latest inferences of primordial el¬ 
ement abundances, 2) using the latest measurements of 
the CMB power spectrum measurements from Planck , 
3) having improved calculations of CMB spectral distor¬ 
tions. 


We also highlight the importance of a region in param¬ 
eter space near m = IMeV and r = 100 ms which we 
term the MeV-ALP window. This region is interesting 
because it previously has evaded all known constraints 
(also noted in [12, 15]) and, as we show, can correspond 
to a particular axion model we will call the DFSZ-EN2. 
Additionally, in this mass window the symmetry breaking 
scale is much lower than the often-considered “invisible” 
axion models. A main conclusion of this paper is that the 
possibility of the DFSZ-EN2 or any other ALP hiding in 
plain sight in the MeV-ALP window is ruled out by the 
newer bounds presented here. We show, however, that 
these bounds are model dependent and the MeV-ALP 
window can be reopened if there is other exotic radiation 
present. 


The paper is structured as follows. In Sec. II we dis¬ 
cuss in more detail the scenario and its implications for 
cosmology. Sec. Ill describes the new and tabulated con¬ 
straints which we use. In Sec. IV A we further discuss 
the MeV-ALP window and in Sec. IV D we give forecasts 
for future probes. 
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II. THE SCENARIO 

We begin by discussing the cosmological impact of 
ALPs, which we define as any particle with a mass 
and two-photo coupling g Following standard conven¬ 
tions in the ALP literature, the effective Lagrangian is 

C = - 9 -f<t>F, v F»\ ( 1 ) 


number of relativistic species, N e fi. As usual, this is taken 
so that in the limit of complete neutrino decoupling prior 
to electron-positron annihilation, each neutrino species 
(with antineutrinos) contributes 1 to N e g, making the 
total relativistic energy density, 
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where F is the electromagnetic field strength tensor, F 
its dual, and (f) the ALP field. We often describe the two 
dimensional parameter space with and, in place of 
<? 0 7 , the lifetime for decay into photons 
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Two processes drive the cosmological evolution of the 
ALP energy density. The first is the Primakoff inter¬ 
action which allows for conversion between photons and 
ALPs in the presence of a charged particle q , via 7 q -f->- <pq. 
Because it is a four-point interaction, the scattering rate 
for the Primakoff process will depend on the density of 
scatterers (in this case charged particles). At early times, 
this density grows faster than the Hubble rate, mean¬ 
ing the Primakoff process will always begin in equilib¬ 
rium and freeze out at later times. The second process is 
the direct two-photon interaction, 77 o (f>. Conversely, 
this three-point interaction has no dependence on scat¬ 
terers, and will always begin out of equilibrium then re¬ 
equilibrate at later times. 

Qualitatively, the ALP scenario depends strongly on 
the time ordering in which ALPs 1) freeze-out from the 
Primakoff interaction, 2) become non-relativistic, and 3) 
recouple via the two-photon interaction. The details of 
how these events are controlled by the two free parame¬ 
ters forbids certain orderings, and in fact the vast ma¬ 
jority of solutions fall into one of just two cases. If 
there is a gap between freeze-out and recoupling dur¬ 
ing which the ALPs become non-relativistic, meaning 
Tf 0 > TO^ > T re , there is an out-of-equilibrium decay. In 
this case, upon becoming non-relativistic ALPs cease to 
track their equilibrium abundance and instead increase 
in energy density relative to the plasma. Decay hap¬ 
pens when the two-photon interaction becomes effective, 
which here is controlled by only the ALP lifetime, inde¬ 
pendent of the mass. On the other hand, if ALPs be¬ 
come non-relativistic only after recoupling, T re > m^, 
they will track their equilibrium abundance throughout 
decay, a scenario which might better be called a “Boltz¬ 
mann suppression.” This suppression occurs when the 
temperature reaches the mass, independent of the life¬ 
time. Although some subtleties can occur if other events 
reheat the plasma while the ALPs are decoupled, to a 
large degree only these two in- and out-of-equilibrium 
scenarios are important. 

In discussing the cosmological impact of ALPs, it is 
useful to define two quantities. The first is the effective 


The other quantity is the baryon-to-photon ratio 77 , or 
equivalently the energy density in baryons flfc/i 2 , 

7 = — = 2.74 x l(T 8 ft 6 /i 2 . (4) 

Tlry 

We sometimes superscript these quantities with BBN or 
CMB depending on the epoch at which they are eval¬ 
uated, although we note there is no exact definition as 
they can change with time in this scenario. 

ALP decays can only reduce or leave unchanged the 
value of A^ e ^ MB . If the decays occur while the neutri¬ 
nos are still fully coupled then A^ t MB does not change. 
Later decays increase the temperature of the plasma rel¬ 
ative to the neutrinos and thus decrease A^ MB . The 
value of iV e BBN can also be decreased in the same man¬ 
ner if decay happens before BBN, although this region 
in parameter space is fairly small. A much larger region 
of parameter space corresponds to ALPs decaying after 
BBN and increasing A”JjP N by simply contributing to the 
relativistic energy density during BBN. ALP decays af¬ 
fect 77 similarly to N e g by causing its value to be less 
after the decay. The difference is that with ? 7 CMB held 
constant, 77 before the decay is now increased. Assuming 
that the physics of BBN is unchanged, both an increase 
in 77 BBN and an increase in -/VjjP N serve to increase the 
amount of primordial helium produced. For deuterium, 
however, there is partial cancellation, leaving the primor¬ 
dial abundance only slightly reduced. This cancellation is 
key to allowing the ALP scenario at all in light of the very 
tight constraints on primordial deuterium (see Sec. Ill B) . 
The increase in primordial helium coupled with the de¬ 
crease in A^ MB moves along a degeneracy direction for 
the CMB constraints, and is hence also generally allowed. 

More quantitatively, we must track the evolution of the 
phase space distribution function /^, which is governed 
by the Boltzmann equation 

^ = (c 9 + c 7 x/; q -/ 0 ), ( 5 ) 

where the /’s are comoving, /^ q is a Bose-Einstein distri¬ 
bution, and a dependence on the comoving momentum p 
is omitted for brevity. When electrons and positrons are 
the only charged particles present, the scattering rate for 
the Primakoff interaction, C q , is given by 


C q 


2 ^ 
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FIG. 1. Key regions and contours in the mass-lifetime pa¬ 
rameter space according to the analytic approximations in 
Sec. II. As in other plots in this paper, dashed lines corre¬ 
spond to the temperature at neutrino decoupling, dot-dashed 
lines the start of BBN, and dotted lines the end of BBN. Blue 
lines show contours of constant Primakoff freeze-out temper¬ 
ature, Tf 0 , and black lines show contours of constant two- 
photon re-equilibration temperature, T re . The line T re = 
divides two regions A and B. Region A is vertically hatched 
and corresponds to out-of-equilibrium decays. Region B is 
cross hatched and corresponds to in-equilibrium decays. Con¬ 
stant decay-time contours in region A are T re = const whereas 
they are = const in region B. Region C has no hatching 
and corresponds to decays before neutrino decoupling, where 
ALPs leave no cosmologically observable traces. The line 
7fo=QCD leaves a sharp feature on cosmological constraints 
as g * changes suddenly during this phase transition. 


where n ep is the number density of electrons and 
positrons, E is the ALP energy, and m 7 = eT /3 is the 
plasmon mass in an electron-positron plasma [16]. As 
shown by Cadamuro and Redondo [11], an accounting of 
other charged particles gives an approximation for the 
Primakoff freeze-out temperature of 


Tf 0 
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GeV, 


( 7 ) 


where is the number of relativistic degrees of freedom 
in the plasma, and g q the number of charged relativistic 
degrees of freedom [17]. 

The two-photon scattering rate, C 7 , can be calculated 
exactly at all times, 


C 1 = 


Et^ 


2T sinh^±E 
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Here we have ignored the plasmon mass, which inhibits 
decays when < 2 m 7 , as this effect is not important in 


the parameter range of interest. We have also assumed 
that photon rethermalization is instantaneous, which is 
an excellent approximation at early times. We consider 
thermalization at later times in Sec. IIIAl. 

In practice, we begin evolving the Boltzmann equa¬ 
tions at a temperature To when electrons and positrons 
are the only remaining charged particles and we can 
thus use Eqn. 6 for the Primakoff scattering rate. If 
the Primakoff process has yet to freeze out at T 0 , as 
given by Eqn. 7, we take take as the initial condi¬ 
tions for a Bose-Einstein distribution at temperature 
To. If freeze-out is earlier, we use the conservation of 
comoving entropy to calculate the increase in photon 
temperature after Primamkoff freeze-out. The result is 
that ALPs are now at a reduced relative temperature of 
To{g*{T 0 ) / g^iTfe)) 1 / 3 . 

We now describe key regions in the ALP parameter 
space using Fig. 1 as a guide. An important quantity 
governing the ALP evolution is the temperature, T re , at 
which the two-photon interaction re-equilibrates. Eqn. 5 
shows that the ALP distribution at a given momen¬ 
tum p becomes equal to its equilibrium value /^ q when 
the scattering rate, C 7 ip ), is on the order of the Hub¬ 
ble constant. We consider only those momenta which 
contribute dominantly to the total energy density, since 
this is the quantity we are interested in. For m,p T, 
these are p « 0. In this case, C 7 vastly simplifies to 
1 /T 0 7 and thus recoupling occurs when the Hubble time 
reaches the lifetime. For < T, the important mo¬ 
menta are instead E ss p ss T. Ignoring the term in 
brackets in Eqn. 8 which can be shown to depend only 
logarithmically on temperature for these momenta, and 
using H = 1.66^TT 2 , we arrive at (in Planck units), 
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Fig. 1 shows constant T re contours in black for three val¬ 
ues of T re corresponding to neutrino decoupling, and the 
start and end of BBN. Here and throughout this paper, 
we take these temperatures to be IMeV, 200 KeV, and 
20 KeV respectively. Fig. 1 also shows constant mass con¬ 
tours at these three values. The line implicitly formed by 
T re = TO 0 divides regions A and B where the ALP be¬ 
comes non-relativistic before and after the two-photon 
interaction re-equilibrates, respectively. Note that in re¬ 
gion A we have rri^ > T re , which is only part of the 
requirement for an out-of-equilibrium decay. The other 
is that Xfo > 777 , 0 , however it turns out that this is al¬ 
ways satisfied for any combination of mass and lifetime 
in region A; that is to say, ALPs can never decay by the 
Primakoff process alone. 

The evolution of the energy densities in the relevant 
components of the plasma for a typical in-equilibrium 
decay (region B) is shown in the third panel of Fig. 2. 
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If ALPs decay after neutrinos are decoupled, conserva¬ 
tion of comoving entropy implies that the temperature 
of the neutrinos relative to the photons after the ALP 
decay and after electrons and positrons have annihilated 
is (4/13) 1 / 3 , as compared to (4/11) 1 / 3 in the standard 
scenario. This means that Ad^ MB will be reduced by 
(13/11) 4 / 3 and, at fixed ? 7 CMB , the value of rj prior to the 
decay is increased by a factor of 13/11. Assuming three 
neutrinos, this gives ~ 2.44. Because the exact 

timing with respect to electron-positron annihilation is 
unimportant for the final temperature, the entirety of re¬ 
gion B shares this same value for Ad^ MB . The BBN data, 
however, are sensitive to the exact time of decay via sen¬ 
sitivity to -/Vjf BN and ?7bbn- Since decay time when in 
equilibrium is controlled only by m^, we find character¬ 
istic constant-mass contours in this region in the BBN 
constraints in Figs. 3, 4, and 5. 

The second panel of Fig. 2 shows instead a typical out- 
of-equilibrium decay (region A). In this region an impor¬ 
tant quantity is the fractional increase in the energy of 
the photons (or equivalently the decrease in N e g) once 
the ALP decays at t ~ r^ 7 . Because the ALP here is 
non-relativistic, this is, 


1 

-^eff 
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where we have assumed that the ALP does not remain 
non-relativistic for too long before decay hence the uni¬ 
verse is radiation dominated and a(t) ~ \Jt. This as¬ 
sumption is true for models right on the edge of the al¬ 
lowed region, thus we find characteristic contours of con¬ 
stant mijjy/ffry for cosmological constraints in region A as 
seen in Figs. 3, 4, and 5. 

In both regions A and B there are characteristic con¬ 
tours arising from phenomena which depend on the frac¬ 
tional energy injection after a certain reaction has frozen 
out, i.e. after a certain temperature. For example, CMB 
and BBN constraints are only sensitive to energy injected 
after neutrino decoupling, as any earlier injection is “in¬ 
visible” because it is rethermalized among all compo¬ 
nents. We will also consider bounds from CMB spectral 
distortions, which depend on the energy injected only af¬ 
ter the freeze-out of reactions which can bring the CMB 
spectrum back into chemical equilibrium. 

In region B, fractional energy injection is independent 
of either mass or lifetime, and the amount after a certain 
temperature is controlled only by the mass. On the non¬ 
equilibrium side the story is slightly more subtle. The 
fractional energy increase at a certain time or scale fac¬ 
tor depends on the mass, as per Eqn. 10. However, we 
cannot simply assume a ~ 1/T because the ALP decay 
alters this relation. In particular, the Friedmann accel¬ 
eration equation shows that it does so in a way which 
exactly cancels the dependence, leaving only sensi¬ 
tivity to r^ 7 . Essentially, more massive ALPs lead to 
more total energy injection, but delay the time it takes 
to reach a certain temperature, leaving the same amount 
of energy injected after that temperature. Thus, contours 


of constant fractional energy injection are m = const in 
region B and r = const in region A. We note that this 
is the same as a contour of constant decay time. The 
final key region in Fig. 1, region C, is delineated by such 
a contour, and corresponds to zero energy injection after 
neutrino decoupling and thus no cosmologically observ¬ 
able imprints. 

The Primakoff freeze-out temperature plays a smaller 
role in the ALP evolution than the recoupling tempera¬ 
ture, although we briefly mention two effects stemming 
from events happening in the gap between freeze-out and 
recoupling (in regions of parameter space where the gap 
exists). We first note that if no reheating of the plasma 
occurs in this gap and ALPs are relativistic, they still 
track their equilibrium abundance even though they are 
decoupled. One possible reheating occurring during the 
gap is the QCD phase transition which imprints a sharp 
feature our constraints. During this transition the tem¬ 
perature of the ALPs roughly doubles along with the rest 
of the plasma if they are still coupled, leading to an en¬ 
ergy density roughly twenty times larger. Another possi¬ 
bility is reheating from electron-positron annihilation, a 
case in which the ALPs can re-equilibrate relativistically 
thereafter. Such a scenario happens near the top-left part 
of Fig. 1 and is shown in the fourth panel of Fig. 2. For 
deuterium and helium constraints shown in Fig. 5, near 
this region we can find both lines of constant T re and Xf 0 . 


III. CONSTRAINTS 
A. Cosmic Microwave Background 

1. Frequency Spectrum 

The measurement of the CMB frequency spectrum by 
COBE/FIRAS places very tight bounds on spectral dis¬ 
tortions away from a black-body spectrum [18], limiting 
possible energy injection into the plasma [19]. 

The effects of energy injection depend crucially on 
when it occurs, with the time-line roughly divided into 
three eras. In the earliest era, reactions that change pho¬ 
ton number are fast and any injection of photon energy 
is quickly rethermalized. This leads to only an adjust¬ 
ment of the temperature, hence this is called the “T era”. 
Decays during this era do not cause any distortions in 
the frequency spectrum observed today and are thus al¬ 
lowed. At about T ~ 750 eV, double Compton (DC) 
and Bremsstrahlung (BR) reactions become too slow to 
significantly change photon number. Injected photons 
now lead to a spectrum with a chemical potential, giving 
this the name the “n era”. Finally, Compton scattering 
freezes out at T « 25 eV, making even the redistribution 
of energy inefficient. Initially we are left with a Comp- 
tonized spectrum described by a y parameter, and at 
even later times photons from decay can be seen directly, 
contributing to the extra-galactic background. We do not 
consider y distortion bounds in this work as they rule out 
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FIG. 2. The evolution of the energy densities in the various 
components of the universe for different scenarios which have 
similar decay time. The temperature of the photons today is 
held fixed and the y-axis units are such that the final value 
of the neutrino line is the value of N^ JLB . As in other plots 
in this paper, dashed lines correspond to the scale factor at 
neutrino decoupling, dot-dashed lines the start of BBN, and 
dotted lines the end of BBN. The dashed red line is not ac¬ 
tually a component, but is shown for illustrative purposes; it 
is the equilibrium ALP energy density (that is, the energy 
density APLs would have if they were in chemical and kinetic 
equilibrium with the photons). As per Eqn. 5, interactions 
serve to always drive the ALP energy density towards equi¬ 
librium. The plots labeled A and B correspond to the same 
regions in Fig. 1. 


only very late decays which are outside of the parameter 
range of interest. 

As shown by Hu and Silk [20], if p(z) -C 1 for the du¬ 
ration of the energy injection and if Compton scattering 
is fast compared to DC and BR, the evolution of p obeys, 


dp dp s /I 1 

dt dt ^ V^dc ^BR 


( 11 ) 


where 1 /£dc and 1 Abr are the scattering rates given by, 


toe _ 1-9 x 10 35 / T \ 9/2 

sec VLbh 2 (l — Yp/2) \ Kelvin/ 

t BR _ 8.9 x 10 26 / T \ ^ 13/4 

sec (Hf,/i 2 ) 3 / 2 (l — Yp/2) \ Kelvin/ 

and dp s /dt is a source term given by, 


dp s 


3tt 4 C(3) 
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( 12 ) 

(13) 


(14) 


which is the chemical potential induced by adding a 
comoving photon energy density dp 7 and number dn 7 . 
Since ALPs decay to two photons, we have dp 7 = — dp$ 
and dn-y = —2dn^,. 

In practice, we first evolve Eqn. 5, which assumes in¬ 
stantaneous thermalization. Then we take the resulting 
solution for dpdn^, and expansion history, and nu¬ 
merically integrate Eqn. 11 to arrive at the actual ther¬ 
malization history. This procedure is correct to first or¬ 
der in the energy injection. We then compare with the 
CORE/FIRAS data which gives a 95% upper bound of, 

\p\ < 9 x 10” 5 (15) 


The excluded region is shown in orange in Fig. 3 and 
is labeled “CMB Spectral Distortions”. Because Eqn. 14 
shows the source term for p depends on the fractional 
energy density injected into the plasma after DC and 
BR freeze-out, the arguments of Sec. II imply that the 
exclusion region should be one of constant decay time, 
and indeed it is parallel to v-dec and BBN start/end 
lines. The exact division between allowed and disal¬ 
lowed scenarios is close to the transition from the T to 
the p era, as expected because almost any chemical po¬ 
tential is ruled out while almost any injection in the T 
era is unobservable. These same constraints were cal¬ 
culated by Cadamuro and Redondo [11] under a further 
approximation of Eqn. 11 where DC and BR are taken 
to be infinitely fast until they instantaneously freeze-out 
at T ss 750 eV. In the “previous constraints” panel of 
Fig. 4 we reproduce their result, showing that qualita¬ 
tively this is a very good approximation. The use of 
Eqn. 11 becomes more important, however, for constrain¬ 
ing scenarios with even smaller chemical potentials gener¬ 
ated deeper into what is currently called the T-era. Two 
future missions which are predicted to reach such sensi¬ 
tivity are PIXIE and COrE [21, 22], for which we give 
forecasts in Sec. IV D. 
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FIG. 3. Exclusion regions in the ALP mass-lifetime parameter space. The dashed and dotted lines labeled “v dec” (neutrino 
decoupling), and “BBN start/end” correspond to particles which decay at these particular times (with decay here arbitrarily 
defined as when maximum energy injection occurs). The two thick dashed lines are the consistency relations for two particular 
axion models (see Sec. IV A). The CMB, D/H, and Y v regions are excluded at 3cr, the Collider and Beam Dump regions are 
excluded at 2a, and the SN1987a and HB Stars regions are less formal, rough bounds (see Sec. IIID). 


2. Angular Power Spectrum 

Measurements of CMB anisotropies have been recently 
improving, both from the ground [23-26] and from space 
[27, 28]. Better angular resolution and lower noise have 
tightened up small-scale constraints where the CMB is 
most sensitive to changes in N e g and Y p , both of which 


are altered by the decay of ALPs. 

A fully general treatment would include ALPs in the 
set of Boltzmann equations for calculating the CMB 
power spectrum, but it turns out this is not necessary 
for the scales that are currently well measured. For these 
scales, all of the physical effects of ALPs are in fact iden¬ 
tical to changes in N e g and Y p , as long as we assume 
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Previous Constraints This Work 



FIG. 4. A comparison of exclusion regions from previous works (left panel) and those presented here (right panel). The right 
panel is identical to Fig. 3. 


adiabatic initial conditions. This is essentially because 
decays must happen early enough, as enforced by the 
spectral distortion bound discussed in the previous sec¬ 
tion, which requires the decay happen by T ss 750 eV 
or equivalently z ~ 3 x 10 6 . The angular scales well 
constrained by CMB measurements, roughly t < 3000, 
correspond to physical scales which do not begin to enter 
the horizon until about z ~ 3 x 10 5 . At this point two 
things are different in the ALP scenario as opposed to the 
standard case 1 ) the amplitude of neutrino density per¬ 
turbations upon horizon entry is reduced relative to the 
photons and 2) the expansion rate is different. However, 
both are exactly captured in the standard scenario by 
changing N e g. No other scale-dependent changes to spa¬ 
tial perturbations are possible because the relevant scales 
are still outside of the horizon by the time of the decay. 
Finally, we note that the altered helium abundance is 
taken as an input to the CMB spectrum calculation, and 
in the ALP scenario it is now just at a different value. 

In their work, Cadamuro and Redondo [11] took as 
a CMB constraint a lower bound on N e g as given by 
WMAP7. The Planck data tighten this constraint and 
are also sensitive to Y p . We use the joint constraint 
on these two parameters given by the combination of 
Planck+ WP+highL from [27], approximating the like¬ 
lihood as Gaussian and taking just the mean and co- 
variance. In practice, we first calculate JV^ MB and Y p 
for a given mass and lifetime (here ignoring uncertain¬ 
ties in Y p due to 77 and nuclear reaction rate uncertain¬ 
ties which are unimportant at the level of the CMB 
constraint). We then calculate the y 2 of these values 
against the mean and covariance of N^g AB -Y p from the 
Planck chain. Masses and lifetimes excluded at > 3er, 
or equivalently y 2 > 11.8, are shown in cyan in Fig. 3 


and labeled “CMB Anisotropies.” This region covers a 
large part of the parameter space corresponding to out- 
of-equilibrium decays. These decays tend to greatly re¬ 
duce N™ B because ALPs can live long enough to be¬ 
come non-relativistic without becoming Boltzmann sup¬ 
pressed, thus greatly increasing their total energy injec¬ 
tion. The majority of the remaining parameter space 
corresponds to in-equilibrium decay, which is allowed be¬ 
cause here there is a lower limit of lV/j/ MB ps 2.44 as 
discussed in Sec. II, compatible with CMB constraints. 
Fig. 6 shows typical values for N^g AB and Y p in the ALP 
scenario as compared to the CMB constraints. 


B. Primordial Abundance Inferences 

The primordial production of the light nuclides dur¬ 
ing BBN occurs from t ~ 1 sec to ~ 3 minutes, corre¬ 
sponding to T ~ 1 MeV to 100 keV [reviewed recently 
in, e.g., 30-32]. Observations of primordial light element 
abundances thus probe new physics at play during this 
epoch [e.g., 33]. For our case of an ALP, light element 
production is affected by changes to the cosmic expan¬ 
sion rate during BBN and to the extrapolation of 77 back 
from the CMB epoch. For the majority of the parameter 
space where ALPs decay after neutrino decoupling, they 
are still present during BBN and hence increase N B ™ 
and the expansion rate. Additionally, because the decay 
decreases rj, fixing 77 cmb to the observed value generally 
leads to an increased ?7bbn- Both effects serve to increase 
Y p , but partially cancel for D/H and 7 Li. The increase in 
ATgP N increases D/H while the increase in 77 has the op¬ 
posite effect, with the latter about twice as large, leaving 
an overall reduced D/H. For 7 Li it is instead the former 
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FIG. 5. The colored contours show the prediction for each of the labeled quantities as a function of different values of ALP 
mass and lifetime. The dotted/dashed/solid lines give 1/2/3 u contours given the measurements for these quantities discussed 
in Sec. III. No lines are visible on the lithium plot because the entire parameter space is excluded at > 3cr (our scenario does not 
alleviate the lithium problem). We do not give contours for the N™ B plot because the CMB constraint is highly degenerate 
with Y p . For the D/H panel, the colored contours are calculated assuming a best-fit r/ from the CMB, and uncertainties in r/ 
and nuclear reaction rates are taken into account in producing the a contours (see Sec. IIIB for discussion). 


which wins out. The light-element trends in the mass¬ 
lifetime planes of Fig. 5 bear out these expectations, as 
we now see in detail. 

In practice we have modified the AlterBBN code of 
Arbey [34] to include changes to the expansion history 
and r] due to ALPs. Our code assumes the photon spec¬ 
trum is instantaneously rethermalized, in effect ignoring 
the possibility that high energy photons from the de¬ 


cay can break apart already-formed nuclei. Bounds due 
to this phenomenon constitute so called “photo-erosion” 
bounds, discussed in e.g. [35] and references therein. We 
will consider them separately at the end of this section. 

Deuterium is observable at z ~ 3 in QSO absorption 
systems, via the ~ 82 km/s isotope shift between D and 
H Lyman absorption lines. Recent D/H measurements 
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FIG. 6 . The contours show the 1- and 2-cr confidence re¬ 
gions for N^ AB and Y p from PAmcfc+WP+highL. The dotted 
lines give the 1 -a constraint on Y p from Aver et al. [29]. The 
dashed line is the relation if standard BBN is assumed, and 
the dot along this line corresponds to the standard value of 
A r ( § f MB = 3.046. Colored points show values of and Y p 

arising from ALP masses and lifetimes sampled from a grid 
over the entire region shown in Fig. 3. They are colored by 
m 0%/ 7 >7 which is an important quantity for the CMB con¬ 
straint since it controls the fractional energy injected into the 
photons (Eqn. 10). The maximum value for in-equilibrium 
decays of N^ JlB « 2.44 gives the sharp cutoff visible above. 
Points along the standard BBN consistency line arise from 
decays happening between neutrino decoupling and the be¬ 
ginning of BBN, and correspond to the island of low helium 
visible in Fig. 5. 


have been reported [36], 

- = (2.53 ± 0.04) x 10 5 , (16) 

H 

which represents a factor ~ 3 improvement in precision. 
These bounds are now so tight as to place the measure¬ 
ment errors on level footing with uncertainties associated 
with nuclear reaction rates and with a determination of 
77 from the CMB. To account for these uncertainties, we 
first consider the joint likelihood for the CMB, D/H, and 
nuclear reaction rate measurements, which can be writ¬ 
ten as 

[DH(m f T 0 7 , 77 , at) - DH]“ 
log A — 2 

zo MEAS 

+ log .£cmb {yh N e g , n) 

+ log^NUCL^i, ^0 (17) 

where DH±ctmeas = (2.53±0.04) x 10 5 as per Eqn. 16, 77 
and N e g are evaluated at the CMB epoch but we omit the 
label for brevity, N e g = N e g (m^, r^ 7 ) is uniquely set by 
the mass and lifetime, a* are parameters describing the 
nuclear reaction rates, and Q and fl 1 are any remaining 
cosmological and nuisance parameters. N e g is important 
and appears explicitly because it is both dependent on 
the ALP parameters and its measurement from the CMB 


is significantly degenerate with 77 . We next analytically 
marginalize over all parameters other than and r ^ 7 
under the assumption that these other parameters have 
Gaussian posterior likelihoods and that D/H depends lin¬ 
early on 77 . This gives 


DH(i ti #, T 0 7 ,77 + ra v °/ f cff , cm) - DH 

— log£ =-- ——2 - 2 -~— 2 —j-~ 

z Fmeas + ct nucl + °etaJ 

(18) 

with 

2 / dDH \ 2 2 n 

o-eta = (! “ r ) ( 19 ) 

where 77 ± a v and N e g ± <T7v eff are the mean and standard 
deviation of the posterior likelihoods from the CMB with 
all other parameters marginalized over, and r is the cor¬ 
relation coefficient between 77 and N e g. The presence of 
r in this equation can be understood by considering the 
r = 1 case, which would imply that CMB measurements 
could turn a fixed N e g into a perfect determination of 77 ; 
the quantity above at which the D/H prediction is evalu¬ 
ated, fj + ra v (N e g — N e g) is the mean of this deter¬ 
mination. Because in our case N e g is fixed by the mass 
and lifetime, it would mean 77 is also fixed, leading to no 
extra uncertainty in D/H. In reality, we find r « 0.4 from 
the Planck measurements. We make one further approx¬ 
imation which is that neither the D/H derivative nor the 
nuclear reaction rate uncertainty depends on the values 
of mass and lifetime or the fact that 77 evolves with time 
in the ALP scenario, which we have checked is sufficient. 
We find ctnucl = 4.5 x 10“' using AlterBBN which takes 
the a.i to be principal components in the nuclear reac¬ 
tion rate parameter space [37]. Numerically evaluating 
the D/H derivative and taking posterior likelihoods from 
Planck+WP+highL, we find cteta = 6.9 x 10“ '. When 
added in quadrature these lead to an effective deuterium 
constraint of 


- = (2.53 ±0.091) x 10 5 (20) 

H 

which is meant to be compared to a theoretical prediction 
calculated for the particular values of 77 and oti given in 
Eqn. 18. 

The effects of these D/H constraints on the (m 0 ,T 0 7 ) 
plane appear in Fig. 5. We see that the effect of an ALP 
is always to decrease D/H due to the ALP’s effective in¬ 
crease of 77 BBN winning out over the increase in 7V^ N . 
Moreover, we see that the high precision of the D/H mea¬ 
surements leads to a tight constraint on the ALP space 
in all regions where the decays occur after neutrino de¬ 
coupling. Indeed, D/H is now a very powerful probe of 
ALPs. 

The primordial 4 He abundance is inferred astronom¬ 
ically from observations of emission spectra of highly 
ionized gas in primitive nearby dwarf galaxies, i.e., in 
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low-metallicity extragalactic HII regions. The primor¬ 
dial abundance is traditionally inferred by extrapolation 
to zero metallicity. To derive helium and metal abun¬ 
dances from the observed spectra requires characteriza¬ 
tion of the thermodynamic properties of the emitting gas 
(i.e., temperature, density). The analysis of [29] derives 
these quantities simultaneously in a self-consistent man¬ 
ner, and finds a primordial abundance 

Y p = 0.2465 ± 0.0097 (21) 

where the uncertainty is quantified with an MCMC anal¬ 
ysis. We adopt this as our fiducial primordial 4 He con¬ 
straint. Using a similar data set, Izotov and Thuan [38] 
give a helium constraint of, Y p = 0.254 ± 0.003 where 
the errors are derived in a less conservative manner. In 
Fig. 5 we see that, as expected, the effect of an ALP is to 
increase Y p for almost all of the parameter space where 
the decays occur after neutrinos decouple. However, the 
constraints are not as strong as those of D/H. There is 
an island of parameter space around > 100 MeV 
and T 0 7 ~ 1 sec where Y p decreases. This region corre¬ 
sponds to decays happening between neutrino decoupling 
and the start of BBN. Here we have ?7bbn = ?7cmb and 
N e g < 3, the latter of which serves to decrease Y p . It 
is interesting to note this low helium region does not ex¬ 
tend along the entire constant decay-time contour, cut¬ 
ting off once we enter the in-equilibrium decay side. This 
occurs because in-equilibrium decays reduce the ALP en¬ 
ergy density more slowly than do out-of-equilibrium ones, 
and thus the entire decay cannot fit in the short time be¬ 
tween neutrino decoupling and BBN. 

Finally, the primordial 7 Li abundance is inferred from 
observations of the atmospheres of low-metallicity (ex¬ 
treme Population II) stars in the Galactic stellar halo. 
Down to a metallicity of (Fe/H) ~ 10 _2 8 (Fe/H)Q, these 
stars have lithium abundances that are the same to 
within a small scatter consistent with observational er¬ 
rors. The independence of Li with Fe in this “Spite 
plateau” indicates that lithium is primordial [39], and 
implies a primordial abundance 

^ = (1.6 ±0.3) x 1CT 10 (22) 

H 

[40]. At lower metallicity, however, the Li/H abundance 
scatter increases dramatically, but always below the Spite 
plateau value. This suggests that in these very metal- 
poor stars some lithium destruction has occurred; the 
reason for this remains unclear. 

The astronomically-inferred lithium abundance in 
Eqn. 22 is inconsistent with the primordial value ex¬ 
pected from standard BBN theory combined with CMB 
determinations of g. The observed Li/H value is low at 
the ~ 5 a level. This is the “lithium problem” [reviewed 
in, e.g., 32]. Stellar astrophysics uncertainties may be 
the origin of the problem, but solutions to date require 
fine tuning and do not explain the observed Li/H “melt¬ 
down” at very low metallicities. A more radical and in¬ 
triguing solution is the presence of new physics during 


or after BBN. The challenge for such scenarios is to re¬ 
duce lithium substantially without drawing other light 
elements-particularly deuterium-from their concordant 
primordial abundances. The ALP scenario tends to ag¬ 
gravate the problem by increasing lithium slightly, as seen 
in Fig. 5, and as expected due to the ALP effect on rj. 
Thus, awaiting a resolution to the lithium problem, we 
do not consider lithium bounds. 

In closing we note that our calculations have neglected 
the effects of photoerosion of the light elements. This oc¬ 
curs 1) when the ALP mass exceeds light-element bind¬ 
ing energies > B ~ 10 MeV, and 2) for decay time 
scales long enough so that the decay photons interact 
with light elements before thermalization. This leads to 
some deuterium destruction via yd — > np , but a net pro¬ 
duction due to e.g., y 4 He —> dd. Thus constraints arise 
from D/H, Y pi and 3 He/D [41-43]. These were recently 
computed for purely electromagnetic decays by Cyburt 
et al. [35], assuming that the decay photons provide a 
negligible contribution to the energy density and thus 
expansion rate. In this case, the constraints are only im¬ 
portant for Tx > 10 4 sec and mx > 10' eV, with X the 
decaying particle. This regime shows a Y p drop due to 
photoerosion and a corresponding D/H increase. These 
trends could potentially bring Y p and D/H predictions 
back into agreement with observations, but would require 
a more detailed calculation. Since the regions of parame¬ 
ter where this can happen are already ruled out by CMB 
observations, we ignore the effects of photo-erosion. 

C. Laboratory 

Laboratory bounds on ALPs come from a variety of 
different experimental setups. At lower masses, roughly 
m < eV, some examples include photo-regeneration ex¬ 
periments (“shining light through walls”), microwave 
cavities, and helioscopes [for a review, see 12-14]. For 
the larger masses considered here, the best constraints 
come from electron-positron colliders and beam dumps. 

The presence of the ALP-photon interaction allows for 
the possibility of single-photon final states at electron- 
positron colliders. Early interpretation of searches for 
these events in terms of constraints on the ALP coupling 
was done by Masso and Toldra [ 8 ]. Both Kleban 
and Rabadan [44] and Mimasu and Sanz [15] have fur¬ 
ther shown the ability of current and future colliders to 
improve these bounds. Here we reproduce the constraint 
from LEP given by Mimasu and Sanz [15] of, 

g^ < 4.5 x HT 4 GeV _1 , (23) 

valid in the entire lifetime range considered here. The ex¬ 
cluded region is labeled “Collider” and shown in magenta 
in Fig. 3. 

Additional constraints come from beam dump exper¬ 
iments, where ALPs would be produced in the beam 
dump, penetrate through shielding, then decay to pho¬ 
tons which can be detected by a downstream detector. 
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We use the constraints from Bjorken et al. [45] which 
find that at 95% confidence, 

m^T^-y > 1.4keVsec. (24) 

This is labeled “Beam Dump” and shown in red in Fig. 3. 


D. Globular Clusters and SN1987A 

ALPs offer a new means for energy loss from stars if 
they can both be produced in stellar interiors and have 
sufficiently weak interaction strengths to subsequently es¬ 
cape. In the case of SN1987A, the energy loss can af¬ 
fect the duration of the neutrino pulse from the hand¬ 
ful of neutrinos which were detected, placing constraints 
on the ALP interaction strength [ 8 , 9]. These bounds 
are reproduced in Fig. 3. We note that they assume 
the ALP has only a two-photon coupling, although con¬ 
straints based on other couplings exist. Energy loss can 
also affect the duration of the red giant phase and of the 
horizontal branch, leading to a different observed ratio of 
such stars in globular clusters. We use the bounds from 
Cadamuro and Redondo [11] based on arguments of Raf- 
felt and Dearborn [46] and Raffelt [47]. These are also 
reproduced in Fig. 3. 

IV. DISCUSSION 


well as standard model fields, transform under the new 
1/(1). The KSVZ model [48, 49] has the standard model 
fermions neutral, whereas in the DFSZ model [50, 51] 
they can carry 17(1) charge. These model dependent 
choices in turn affect the axion’s effective photon cou¬ 
pling which arises from fermion loops, ultimately leading 
to a consistency relation between axion mass and photon 
lifetime which can be written as, 



2(4 + z) 
3(1 + z) 


-l 


4~Z 16?r 3 / 2 
1 + z a 


fitWl-n 


-5 


(26) 


with the model dependence captured by the E/N factor. 
In evaluating this relation, we will adopt fixed values of 
m n = 135 MeV, /' n = 92MeV, and 2 = 0.56, ignoring 
small uncertainties that lead to roughly a 10 % uncer¬ 
tainty in the axion mass [10, 52]. 

The KSVZ model has E/N = 0, so that we have, 

= 6.20 x 10 4 (—^f) 5 (27) 

V sec J V eV / 

The DFSZ model we consider here has E/N = 2, hence 
we refer to it as the “DFSZ-EN2” model, with consis¬ 
tency relation, 

= 2.66 x 10 5 5 (28) 

V sec / V eV / 


A. The MeV-ALP Window 


An interesting feature of the exclusion regions prior to 
this work is the allowed window bounded on all sides near 
m ~ 1 MeV and r ~ 100 ms corresponding to an ALP 
decay during BBN. We will henceforth call this the MeV- 
ALP window. It can be seen in the left panel of Fig. 4 
as well as in Hewett et al. [12], Mimasu and Sanz [15]. 
Further interest is driven by the fact that a particle in this 
window could actually be a DFSZ axion, in which case its 
symmetry breaking scale is close to the electroweak scale. 
One of the main conclusions of this work is to show that 
this region is now, in fact, ruled out by the combination 
of CMB+D/H measurements. 

We first briefly review two relevant generic axion mod¬ 
els, referred to as the KSVZ and DFSZ models. Both 
models introduce a new global 17(1) symmetry which is 
approximately broken at some energy scale giving rise 
to an axion with mass m^. The symmetry breaking scale 
is related to the axion mass by non-perturbative effects 
and given by 




y/z TOtt/tt 
1 + Z U 


(25) 


where m n is the pion mass, its decay constant, 2 = 
in u /md the ratio of up to down quark masses. Axion 
models differ in what other new fields are introduced 
to implement the symmetry breaking and how these, as 


These two consistency relations are shown as the dashed 
lines in Fig. 3. Other values for E/N are possible, but 
the DFSZ-EN2 has the distinction of having a partic¬ 
ularly weak coupling because E/N « 2(4 + z)/2>(\ + z) 
and so these terms nearly cancel in Eqn. 26 [53, 54]. This 
weaker coupling means that the DFSZ-EN2 is consistent 
with the collider bounds over a larger range of masses 
as compared to the KSVZ. Ultimately it is that the con¬ 
sistency relation passes through the MeV-ALP window 
which motivates our interest in this model. The lower 
mass limit for the DFSZ-EN2 in the MeV-ALP window 
is around m$ ~ 200 keV, corresponding to //, ~ 30 GeV, 
less than an order of magnitude from the electroweak 
scale u W eak ~ 246 GeV where the axion was initially 
thought to lie. 

While it is interesting that this mass range for the 
DFSZ-EN2 was previously allowed, this part of param¬ 
eter space for the DFSZ-EN2, and more generally the 
entire MeV-ALP window, is now ruled out by the com¬ 
bination of CMB+D/H data. This region corresponds 
to in-equilibrium decays hence it gives A^ MB = 2.44. 
The decay happens essentially in the middle of BBN, in¬ 
creasing A(^ BN and ?;bbn , which, as discussed previously, 
increases Y p and decreases D/H. The decrease in A^ MB 
and increase in Y p moves along the degeneracy direction 
for CMB measurements, and is allowed even by our up¬ 
dated CMB constraints coming from Planck (see Fig. 6 ). 
It is in combination with the D/H constraints that the 
MeV-ALP window is closed, with the best fitting model 
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DFSZ-EN2 case 



FIG. 7. Parameter constraints in the MeV-ALP region of 
parameter space when we also allow extra radiation present 
besides neutrinos and the ALP. Our likelihood includes all 
of the constraints shown in Fig. 3. In the ALP case (top 
panel) the lifetime is marginalized over whereas in the DFSZ- 
EN2 case (bottom panel) it is fixed by the consistency re¬ 
lation (Eqn. 28). The vertical dashed line is a forecast for 
SUPER-KEKB, showing that it could close the remaining al¬ 
lowed parameter window or detect a particle there. We show 
A^ e ff evaluated both prior to BBN when neutrinos, extra ra¬ 
diation, and the ~MeV ALP (which here adds 4/7 to JV e ff) 
contribute, and at the CMB epoch after the ALP has decayed. 
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within the window ruled out at about 3.5<r. If we replace 
the CMB constraint from Planck with previous measure¬ 
ments from the combination of WMAP, ACT, and SPT, 
the window is ruled out at a similar significance. This 
is despite the 20% tighter rj constraint from Planck be¬ 
cause the central value also shifts lower, increasing D/H 
back towards the measured value. Conversely, replacing 
the D/H measurement with previous bounds does open 
the MeV-ALP window again, as seen in the left panel of 
Fig. 4. It is thus the new bounds from Cooke et al. [36] 
that are the key improvement. 


B. A Loophole in the Presence of Extra Radiation 


FIG. 8. (Top) The contours show the 1- and 2-cr confidence 
regions for N^ AB and Y p from Plancfc+WP+highL. The dot¬ 
ted lines give the 1 -a constraint on Y v from Aver et al. [29]. 
The dashed line is the relation if standard BBN is assumed, 
and the dot along this line corresponds to the standard value 
of = 3.046. Colored points show values of N^g lB 

and Y p arising from ALP masses and lifetimes taken from 
the ACDM+AAeff+ALP chain described in Sec. IV B. They 
are colored by m^, which controls decay time and can alter 
BBN but otherwise does not affect N^ AB . (Bottom) Same 
as the top panel, but with 10 5 D/H shown on the x-axis. The 
vertical dotted lines give the 1 -<t constraint from Cooke et al. 
[36]. Points are colored instead by N^g eBBN . Sufficiently tight 
constraints around the standard value (black dot) could rule 
out the ALP scenario even in the presence of extra radiation. 


While the MeV-ALP region is now excluded by the 
CMB+D/H measurements, this result depends on the 
assumption of having no extra radiation besides neutri¬ 
nos and the ALP. In some scenarios, for example as pre¬ 
dicted by the string axiverse [7], it is natural to have 
many ALPs, some of which could also contribute to N e g 
but be light enough to remain otherwise invisible. Moti¬ 
vated by this possibility, we explore constraints when in 
addition to the ALP mass and lifetime, we also allow an 
extra arbitrary addition to N e g. 

The MeV-ALP region is ruled out largely because it 
predicts too low an abundance of primordial deuterium. 
An addition to fV^ BN increases D/H and can bring it 
back into agreement with measurements. The penalty 
is a further increase in Y p , but because the helium con¬ 


straints are not as tight as D/H, an allowed window now 
opens up again. 

We explore this window with Markov Chain Monte 
Carlo (MCMC) [55]. We run two MCMC chains, one 
for the ALP case where both mass and lifetime are free 
parameters, and another for the DFSZ-EN2 case with 
only the mass free and the lifetime given by Eqn. 28. In 
both cases we also leave free the quantity we call A N e g 
which controls any extra relativistic energy density at 
some early time before BBN, meaning that ./Vfg CBBN = 
3 + 4/7 + A N e ff. The ALP contributes 4/7 because it is 
one bosonic degree of freedom and is fully thermalized in 
all regions of parameter space explored by the chain. The 
likelihood includes all of the bounds in Fig. 3. In either 
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the ALP of DFSZ-EN2 cases, we find that the MeV-ALP 
window is again allowed and the best-fitting model con¬ 
sistent with all of the data at < la. 

The mass posterior distributions for the ALP and 
DFSZ-EN2 chains are shown in Fig. 7. In both cases 
masses below 200 keV are excluded. Masses above 1 MeV 
in DFSZ-EN2 case are excluded by the collider bound, 
but are allowed in the ALP case because models can 
evade this constraint by having a smaller photon cou¬ 
pling. We also show a forecast for a next generation 
electron-positron collider SUPER-KEKB after two years 
of integration (discussed in Sec. IV D) which can probe 
down to almost exactly the 200 keV minimum. 

The corresponding likelihoods for N e s are given in the 
right panel. These show that the data accommodate the 
ALP scenario by initially having 7V ( f f / oBBN ~ 4.7 and di¬ 
luting this down to N™ B ~ 3.4 via the ALP decay. Con¬ 
straints on the extra radiation are AN e s = 1.13 ± 0.30. 
Thus, one or (marginally) two extra neutrino-like parti¬ 
cles allow for an ALP in the MeV-ALP window. Alterna¬ 
tively, the decay of such an ALP can hide the existence 
of one or two additional neutrino-like particles from the 
tight CMB constraints which would otherwise rule them 
out. A similar loophole allowing for extra radiation has 
been proposed by Ho and Scherrer [56]. 

We next test the extent to which the data prefer 
these extended models. We perform a simple test us¬ 
ing best-fit x 2 values given in Tab. II. If the \ 2 for the 
extended model decreases significantly as compared to 
ACDM, then roughly that model is preferred. Although 
all bounds from Fig. 3 are included in the fit, we only give 
X' 2 for those which are not hard cutoffs. When using the 
combination of Planck+T)/K+Y p , we find the baseline 
data choice slightly disfavors both the DFSZ-EN2 and 
ALP models. The only case where there is a preference 
for the extend model is in the ALP case when using the 
helium constraint from Izotov and Thuan [38]. Here we 
find an improvement in \ 2 °f 5-81 when we have added 
3 new free parameters, something we expect to happen 
by chance only 12% of the time. If the high helium value 
inferred by Izotov and Thuan [38] is confirmed, then this 
scenario is a natural explanation as it can increase he¬ 
lium compared to the standard value while keeping the 
deuterium abundance and IV^ MB roughly unchanged. 


C. A Simple Expression for Exclusion Bounds 

Given the improved constraints from CMB+D/H mea¬ 
surements, we suggest a simple expression for ALP 
bounds which can be adopted by those who prefer a sim¬ 
pler picture than the many probes shown in Fig. 3. The 
CMB+D/H data alone now essentially rule out any en¬ 
ergy injection after neutrino decoupling, giving allowed 
parameters of 


m ^ 
eV 


> 10 7 and ^ 


< io- 


(29) 


This assumes no extra radiation besides ALPs, and is 
valid roughly until masses become small enough or life¬ 
times long enough that decays happen after CMB last 
scattering. These late decays are analyzed in more detail 
by Cadamuro and Redondo [11], who find approximately 

KM < io 1 or T -* > 10 24 , (30) 

eV sec 

are once again allowed. 


D. Forecasts 

Measurements relevant for placing bounds on ALP pa¬ 
rameters have been recently improving and will continue 
to do so in the near future. It is expected that several 
probes will soon have the sensitivity to further test the 
MeV-ALP window. In this section we compute forecasts 
for some of them. 

Currently CMB anisotropies alone are not enough to 
rule out the MeV-ALP window where there is a maxi¬ 
mum of N e g = 2.44 and an increase in Y p , but which 
lie along the CMB degeneracy direction and are thus al¬ 
lowed. Abazajian et al. [57] show that a Stage-IV CMB 
experiment could measure I\dj/ MB to within 0.02 at 1-a. 
Given such tight constraints, the arguments of Sec. Ill A 2 
may need to be revisited; while it is true that the ALP 
decays before any modes relevant for the CMB enter the 
horizon, the difference is only an order of magnitude in 
scale factor. Assuming any such corrections do not pro¬ 
vide loopholes, if a Stage-IV CMB measurement found a 
value of N*g consistent with the standard value of 3.046, 
the MeV-ALP window would be strongly ruled out [58]. 
For the CMB constraints, however, there will always be 
the possibility of extra radiation exactly canceling the 
dilution due to the ALP decay (as in Sec. IVB). 

Helium and deuterium measurements will continue 
to improve as more systems are discovered and sys¬ 
tematic errors are better understood. Additionally, 
D/H measurements can be significantly improved by 
better measuring nuclear reaction rates in the labo¬ 
ratory. Fig. 8 shows that in the D/H -Y p plane, the 
ACDM+AA/fj+ALP scenario is not continuous with the 
standard model. This allows for sufficiently tight Y p and 
D/H constraints around the standard values to rule out 
the presence of ALPs, independent of assumptions about 
extra radiation. We find that the minimum requirement 
for all points in the ACDM+A V e ff+ALP chain to be 
ruled out at > 3cr by the combination of Y p and D/H 
measurements is a factor of two improvement in the D/H 
error and a factor of three improvement in the Y p error 
bar. We note the former is possible by eliminating the 
uncertainty due to nuclear reaction rates alone. 

On the laboratory side, in Sec. Ill C we used the con¬ 
straint from LEP which limited c^ 7 < 4.5 x 10 -4 GeV -1 . 
If a search for single-photon events were performed us¬ 
ing the entire 1000 fb~ 4 of currently existing KERB data 
and the standard model background was found, forecasts 
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TABLE I. Best-fit parameters 



m [keV] 

t [ms] A^ eBBN 

at-CMB 

10 5 D/H 

Y P 

10 10 Li 7 

ACDM 


3 

3.046 

2.56 

0.247 

4.58 

ACDM+AALff+ALP 

1936 

4.6 4.73 

3.52 

2.46 

0.255 

5.02 

ACDM+AA eff +DFSZ-EN2 

734 

6.2 4.61 

3.30 

2.45 

0.258 

5.15 


TABLE II. Best-fit * 2 



Planck(2) Cooke(l) Aver(l) Izotov(l) Planck+Cooke-|-Aver(4) Planck+Cooke+Izotov(4) 

ACDM 

0.96 

0.10 

0.00 

5.62 

1.06 

6.68 

ACDM+AAeff+ALP 

0.22 

0.57 

0.74 

0.08 

1.53 

0.87 

ACDM+AA eff +DFSZ-EN2 

0.02 

0.91 

1.43 

1.86 

2.35 

2.79 


from Kleban and Rabadan [44] show that the constraints 
could improve to g < 10 _6 GeV _1 . Similar improve¬ 
ment could come from reinterpreting the constraints on 
dark photons from 500 fb _1 of BABAR data given in [59] 
in terms of ALPs. Attempting either of these is outside 
of the scope of this paper, but could make significant im¬ 
provements in the mass bounds shown in Fig. 7. SUPER- 
KEKB, an ongoing upgrade to KERB, plans to improve 
on the integrated luminosity of KERB by a factor of ten 
with two years of integration, and by a factor of fifty with 
ten years [60]. Taking the constraint on g to scale with 
the square root of the integrated luminosity, we find that 
within the first two years, SUPER-RERB can rule out 
the last remaining part of the MeV-ALP window through 
which the DFSZ-EN2 passes (the mass limit forecast is 
shown in Fig. 7). The full ten year forecast is exactly 
enough to close the MeV-ALP window entirely, bring¬ 
ing the collider bound up to g <^ 7 < 10~ 6 GeV -1 where 
the SN1987a constraint begins. This simple forecast is in 
broad agreement with a more sophisticated calculation 
given by Mimasu and Sanz [15]. 

Finally, PIXIE is a proposed mission which would 
greatly improve constraints on CMB spectral distortions 
[21]. Expected bounds on the /i parameter are, 


|/i| < 5 x 1CT 8 . (31) 


V. CONCLUSION 

We have shown how cosmological, astrophysical, and 
laboratory bounds can provide complementary con¬ 
straints in the mass-lifetime parameter space of axions 
and ALPs. We have updated the work of Cadamuro and 
Redondo [11] with constraints from the Planck satellite 
and the latest inferences of primordial D/H and helium 
abundances, and provided a more detailed calculation of 
spectral distortions. The most important change is that 
CMB+D/H constraints now rule out the entire region 
corresponding to decays happening after neutrino decou¬ 
pling but before CMB last scattering. This includes clos¬ 
ing the MeV-ALP region of parameter space which we 
have also shown can correspond to a type of DFSZ axion. 
The presence of additional radiation can relax the exclu¬ 
sion regions and once again allow the MeV-ALP window. 
Although it is allowed in this case, including such a par¬ 
ticle slightly degrades the overall fit to the CMB+BBN 
data if our most robust data combination is used. Alter¬ 
natively, this model can provide a natural explanation for 
a high value of helium such as found by Izotov and Thuan 
[38]. Forecasts for future primordial abundance measure¬ 
ments and for SUPER-RERB are promising; both have 
the ability to test the MeV-ALP window even in the pres¬ 
ence of extra radiation. A detection by either would be 
very exciting. Even the null result, however, would sig¬ 
nify a new level of precision in our understanding of the 
contents of the primordial plasma. 


Bounds due to spectral distortion constitute constant 
decay-time boundaries (Sec. II), and we find that the 
PIXIE forecast moves up the exclusion region by only a 
factor of five in decay time as compared to FIR AS. This 
is much less than needed to reach the MeV-ALP window. 
Evidently it is very difficult to constrain decays happen¬ 
ing very much into the T era using spectral distortions. 


ACKNOWLEDGMENTS 

We thank Xevin Cleary, Francis-Yan Cyr-Racine, 
Rouen Essign, Brent Follin, Nemanja Raloper, Yury 
Rolomensky, Markus Luty, Javier Redondo, and Mc- 
Cullen Sandora for helpful conversations and input. We 
acknowledge support from NSF awards No. 1213801. 


[1] R. D. Peccei and H. R. Quinn, Physical Review Letters 
38, 1440 (1977). 


[2] R. D. Peccei and H. R. Quinn, Physical Review D 16, 
1791 (1977). 











15 


[3] F. Wilczek, Physical Review Letters 40, 279 (1978). 

[4] S. Weinberg, Physical Review Letters 40, 223 (1978). 

[5] Y. Chikashige, R. N. Mohapatra, and R. D. Peccei, 
Physics Letters B 98, 265 (1981). 

[6] G. B. Gelmini and M. Roncadelli, Physics Letters B 99, 
411 (1981). 

[7] A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, 
and J. March-Russell, arXiv:0905.4720 [astro-ph, 
physics:gr-qc, physics:hep-ph, physics:hep-th] (2009), 
arXiv: 0905.4720. 

[8] E. Masso and R. Toldra, Physical Review D 52, 1755 
(1995). 

[9] E. Masso and R. Toldra, Physical Review D 55, 7967 
(1997), arXiv: hep-ph/9702275. 

[10] D. Cadamuro, S. Hannestad, G. Raffelt, and J. Redondo, 
Journal of Cosmology and Astroparticle Physics 2011, 
003 (2011). 

[11] D. Cadamuro and J. Redondo, JCAP 2, 32 (2012). 

[12] J. L. Hewett, H. Weerts, R. Brock, J. N. Butler, B. C. K. 
Casey, J. Collar, A. de Gouvea, R. Essig, Y. Gross- 
man, W. Haxton, J. A. Jaros, C. K. Jung, Z. T. 
Lu, K. Pitts, Z. Ligeti, J. R. Patterson, M. Ramsey- 
Musolf, J. L. Ritchie, A. Roodman, K. Scholberg, 

C. E. M. Wagner, G. P. Zeller, S. Aefsky, A. Afana- 
sev, K. Agashe, C. Albright, J. Alonso, C. Ankenbrandt, 
M. Aoki, C. A. Arguelles, N. Arkani-Hamed, J. R. Ar- 
mendariz, C. Armendariz-Picon, E. A. Diaz, J. Asaadi, 

D. M. Asner, K. S. Babu, K. Bailey, O. Baker, B. Bal- 
antekin, B. Bailer, M. Bass, B. Batell, J. Beacham, 

J. Behr, N. Berger, M. Bergevin, E. Berman, R. Bern¬ 
stein, A. J. Bevan, M. Bishai, M. Blanke, S. Blessing, 

A. Blondel, T. Blum, G. Bock, A. Bodek, G. Bonvicini, 

F. Bossi, J. Boyce, R. Breedon, M. Breidenbach, S. J. 
Brice, R. A. Briere, S. Brodsky, C. Bromberg, A. Bross, 
T. E. Browder, D. A. Bryman, M. Buckley, R. Burn- 
stein, E. Caden, P. Campana, R. Carlini, G. Carosi, 
C. Castromonte, R. Cenci, I. Chakaberia, M. C. Chen, 

C. H. Cheng, B. Choudhary, N. H. Christ, E. Chris¬ 
tensen, M. E. Christy, T. E. Chupp, E. Church, D. B. 
Cline, T. E. Coan, P. Coloma, J. Comfort, L. Coney, 

J. Cooper, R. J. Cooper, R. Cowan, D. F. Cowen, 

D. Cronin-Hennessy, A. Datta, G. S. Davies, M. De- 
marteau, D. P. DeMille, A. Denig, R. Dermisek, A. Desh- 
pande, M. S. Dewey, R. Dharmapalan, J. Dhooghe, 
M. R. Dietrich, M. Diwan, Z. Djurcic, S. Dobbs, M. Du- 
raisamy, B. Dutta, H. Duyang, D. A. Dwyer, M. Eads, 

B. Echenard, S. R. Elliott, C. Escobar, J. Fajans, S. Fa- 
rooq, C. Faroughy, J. E. Fast, B. Feinberg, J. Felde, 

G. Feldman, P. Fierlinger, P. F. Perez, B. Filippone, 
P. Fisher, B. T. Flemming, K. T. Flood, R. Forty, 
M. J. Frank, A. Freyberger, A. Friedland, R. Gandhi, 

K. S. Ganezer, A. Garcia, F. G. Garcia, S. Gardner, 

L. Garrison, A. Gasparian, S. Geer, V. M. Gehman, 
T. Gershon, M. Gilchriese, C. Ginsberg, I. Gogoladze, 

M. Gonderinger, M. Goodman, H. Gould, M. Graham, 
P. W. Graham, R. Gran, J. Grange, G. Gratta, J. P. 
Green, H. Greenlee, R. C. Group, E. Guardincerri, 
V. Gudkov, R. Guenette, A. Haas, A. Hahn, T. Han, 
T. Handler, J. C. Hardy, R. Harnik, D. A. Harris, F. A. 
Harris, P. G. Harris, J. Hartnett, B. He, B. R. Heckel, 
K. M. Heeger, S. Henderson, D. Hertzog, R. Hill, E. A. 
Hinds, D. G. Hitlin, R. J. Holt, N. Holtkamp, G. Horton- 
Smitli, P. Huber, W. Huelsnitz, J. Imber, I. Irastorza, 
J. Jaeckel, I. Jaegle, C. James, A. Jawahery, D. Jensen, 


C. P. Jessop, B. Jones, H. Jostlein, T. Junk, A. L. Kagan, 
M. Kalita, Y. Kamyshkov, D. M. Kaplan, G. Karagiorgi, 
A. Karle, T. Katori, B. Kayser, R. Kephart, S. Kettell, 
Y. K. Kim, M. Kirby, K. Kirch, J. Klein, J. Kneller, 

A. Kobach, M. Kohl, J. Kopp, M. Kordosky, W. Ko- 
rsch, I. Kourbanis, A. D. Krisch, P. Krizan, A. S. Kro- 
nfeld, S. Kulkarni, K. S. Kumar, Y. Kuno, T. Kutter, 
T. Lachenmaier, M. Lamm, J. Lancaster, M. Lancaster, 

C. Lane, K. Lang, P. Langacker, S. Lazarevic, T. Le, 
K. Lee, K. T. Lesko, Y. Li, M. Lindgren, A. Lind¬ 
ner, J. Link, D. Lissauer, L. S. Littenberg, B. Lit¬ 
tlejohn, C. Y. Liu, W. Loinaz, W. Lorenzon, W. C. 
Louis, J. Lozier, L. Ludovici, L. Lueking, C. Lunardini, 

D. B. MacFarlane, P. A. N. Machado, P. B. Mackenzie, 
J. Maloney, W. J. Marciano, W. Marsh, M. Marshak, 
J. W. Martin, C. Mauger, K. S. McFarland, C. McGrew, 
G. McLaughlin, D. McKeen, R. McKeown, B. T. Mead¬ 
ows, R. Mehdiyev, D. Melconian, H. Merkel, M. Messier, 
J. P. Miller, G. Mills, U. K. Minamisono, S. R. Mishra, 

I. Mocioiu, S. M. Sher, R. N. Mohapatra, B. Mon- 
real, C. D. Moore, J. G. Morfin, J. Mousseau, L. A. 
Moustakas, G. Mueller, P. Mueller, M. Muether, H. P. 
Mumm, C. Munger, H. Murayama, P. Nath, O. Naviliat- 
Cuncin, J. K. Nelson, D. Neuffer, J. S. Nico, A. Norman, 
D. Nygren, Y. Obayashi, T. P. O’Connor, Y. Okada, 

J. Olsen, L. Orozco, J. L. Orrell, J. Osta, B. Pahlka, 
J. Paley, V. Papadimitriou, M. Papucci, S. Parke, R. H. 
Parker, Z. Parsa, K. Partyka, A. Patch, J. C. Pati, 
R. B. Patterson, Z. Pavlovic, G. Paz, G. N. Perdue, 
D. Perevalov, G. Perez, R. Petti, W. Pettus, A. Piepke, 
M. Pivovaroff, R. Plunkett, C. C. Polly, M. Pospelov, 
R. Povey, A. Prakesh, M. V. Purohit, S. Raby, J. L. 
Raaf, R. Rajendran, S. Rajendran, G. Rameika, R. Ram¬ 
sey, A. Rashed, B. N. Ratcliff, B. Rebel, J. Redondo, 
P. Reimer, D. Reitzner, F. Ringer, A. Ringwald, S. R.ior- 
dan, B. L. Roberts, D. A. Roberts, R. Robertson, F. Ro- 
bicheaux, M. Rominsky, R. Roser, J. L. Rosner, C. Rott, 
P. Rubin, N. Saito, M. Sanchez, S. Sarkar, H. Schellman, 

B. Schmidt, M. Schmitt, D. W. Schmitz, J. Schneps, 
A. Schopper, P. Schuster, A. J. Schwartz, M. Schwarz, 
J. Seeman, Y. K. Semertzidis, K. K. Seth, Q. Shall, 
P. Shanahan, R. Sharma, S. R. Sharpe, M. Shiozawa, 
V. Shiltsev, K. Sigurdson, P. Sikivie, J. Singh, D. Sivers, 
T. Skwarnicki, N. Smith, J. Sobczyk, H. Sobel, M. Soder- 
berg, Y. H. Song, A. Soni, P. Souder, A. Sousa, J. Spitz, 

M. Stancari, G. C. Stavenga, J. H. Steffen, S. Stepanyan, 

D. Stoeckinger, S. Stone, J. Strait, M. Strassler, I. A. Su- 
lai, R. Sundrum, R. Svoboda, B. Szczerbinska, A. Szelc, 
T. Takeuchi, P. Tanedo, S. Taneja, J. Tang, D. B. Tanner, 
R. Tayloe, I. Taylor, J. Thomas, C. Thorn, X. Tian, B. G. 
Tice, M. Tobar, N. Tolich, N. Toro, I. S. Towner, Y. Tsai, 
R. Tschirhart, C. D. Tunnell, M. Tzanov, A. Upadhye, 
J. Urheim, S. Vahsen, A. Vainshtein, E. Valencia, R. G. 
Van de Water, R. S. Van de Water, M. Velasco, J. Vo¬ 
gel, P. Vogel, W. Vogelsang, Y. W. Wah, D. Walker, 

N. Weiner, A. Weltman, R. Wendell, W. Wester, M. Wet- 
stein, C. White, L. Whitehead, J. Whitmore, E. Wid- 
mann, G. Wiedemann, J. Wilkerson, G. Wilkinson, 
P. Wilson, R. J. Wilson, W. Winter, M. B. Wise, 
J. Wodin, S. Wojcicki, B. Wojtsekhowski, T. Wongjirad, 

E. Worcester, J. Wurtele, T. Xin, J. Xu, T. Yamanaka, 
Y. Yamazaki, I. Yavin, J. Yeck, M. Yeh, M. Yokoyama, 
J. Yoo, A. Young, E. Zimmerman, K. Zioutas, M. Zis- 
man, J. Zupan, and R. Zwaska, arXiv: 1205.2671 [hep-ex, 


16 


physics:hep-ph] (2012), arXiv: 1205.2671. 

[13] R. Essig, J. A. Jaros, W. Wester, P. Hansson Adrian, 

S. Andreas, T. Averett, O. Baker, B. Batell, 

M. Battaglieri, J. Beacham, T. Beranek, J. D. Bjorken, 
F. Bossi, J. R. Boyce, G. D. Cates, A. Celentano, 

A. S. Chou, R. Cowan, F. Curciarello, H. Davoudi- 
asl, P. deNiverville, R. De Vita, A. Denig, R. Dharma- 
palan, B. Dongwi, B. Dobrich, B. Echenard, D. Es- 
priu, S. Fegan, P. Fisher, G. B. Franklin, A. Gaspar- 
ian, Y. Gershtein, M. Graham, P. W. Graham, A. Haas, 

A. Hatzikoutelis, M. Holtrop, I. Irastorza, E. Izaguirre, 

J. Jaeckel, Y. Kahn, N. Kalantarians, M. Kohl, G. Krn- 
jaic, V. Kubarovsky, H. Lee, A. Lindner, A. Lobanov, 
W. J. Marciano, D. J. E. Marsh, T. Maruyama, D. Mc- 
Keen, H. Merkel, K. Moffeit, P. Monaghan, G. Mueller, 

T. K. Nelson, G. R. Neil, M. Oriunno, Z. Pavlovic, S. K. 
Phillips, M. J. PivovarofT, R. Poltis, M. Pospelov, S. Ra- 
jendran, J. Redondo, A. Ringwald, A. Ritz, J. Ruz, 

K. Saenboonruang, P. Schuster, M. Shinn, T. R. Slatyer, 
J. H. Steffen, S. Stepanyan, D. B. Tanner, J. Thaler, 
M. E. Tobar, N. Toro, A. Upadye, R. Van de Water, 

B. Vlahovic, J. K. Vogel, D. Walker, A. Welt man, B. Wo- 
jtsekliowski, S. Zhang, and K. Zioutas, ArXiv e-prints 
(2013). 

[14] K. A. Olive and P. D. Group, Chinese Physics C 38, 
090001 (2014). 

[15] K. Mimasu and V. Sanz, ArXiv e-prints 1409, 4792 
(2014). 

[16] M. Bolz, A. Brandenburg, and W. Buchmiiller, Nuclear 
Physics B 606, 518 (2001). 

[17] The degrees of freedom g q are counted identically to g , 
but weighted by the square of the charge of each species, 
in units of the elementary charge. Eqn. 6 ignores the de¬ 
pendence of the plasmon mass on g q , leading to changes 
in freeze-out temperature as large as ~ 25% [see 11]. 
Since the details of Primakoff freeze-out are not impor¬ 
tant to our main conclusions, this is sufficient. 

[18] D. J. Fixsen, E. S. Cheng, J. M. Gales, J. C. Mather, 

R. A. Shafer, and E. L. Wright, Astrophys. J. 473, 576 
(1996). 

[19] E. L. Wright, J. C. Mather, D. J. Fixsen, A. Kogut, R. A. 
Shafer, C. L. Bennett, N. W. Boggess, E. S. Cheng, R. F. 
Silverberg, G. F. Smoot, and R. Weiss, The Astrophys- 
ical Journal 420, 450 (1994). 

[20] W. Hu and J. Silk, Phys. Rev. D 48, 485 (1993). 

[21] A. Kogut, D. J. Fixsen, D. T. Chuss, J. Dotson, E. Dwek, 
M. Halpern, G. F. Hinshaw, S. M. Meyer, S. H. Moseley, 
M. D. Seiffert, D. N. Spergel, and E. J. Wollack, Journal 
of Cosmology and Astro-Particle Physics 07, 025 (2011). 

[22] T. C. Collaboration, C. Armitage-Caplan, M. Avillez, 

D. Barbosa, A. Banday, N. Bartolo, R. Battye, 

J. Bernard, P. de Bcrnardis, S. Basak, M. Bersanelli, 

P. Bielewicz, A. Bonaldi, M. Bucher, F. Bouchet, 

F. Boulanger, C. Burigana, P. Camus, A. Challinor, 

S. Chongchitnan, D. Clements, S. Colafrancesco, J. De- 
labrouille, M. De Petris, G. De Zotti, C. Dickinson, 
J. Dunkley, T. Ensslin, J. Fergusson, P. Ferreira, K. Fer- 
riere, F. Finelli, S. Galli, J. Garcia-Bellido, C. Gau¬ 
thier, M. Haverkorn, M. Hindmarsh, A. Jaffe, M. Kunz, 
J. Lesgourgues, A. Liddle, M. Liguori, M. Lopez- 
Caniego, B. Maffei, P. Marchegiani, E. Martinez- 
Gonzalez, S. Masi, P. Mauskopf, S. Matarrese, A. Mel- 
chiorri, P. Mukherjee, F. Nati, P. Natoli, M. Negrello, 

L. Pagano, D. Paoletti, T. Peacocke, H. Peiris, L. Per- 


roto, F. Piacentini, M. Piat, L. Piccirillo, G. Pisano, 

N. Ponthieu, C. Rath, S. Ricciardi, J. Rubino Martin, 

M. Salatino, P. Shellard, R. Stompor, L. T. J. Urrestilla, 

B. Van Tent, L. Verde, B. Wandelt, and S. Withington, 
ArXiv e-prints 1102, 2181 (2011). 

[23] R. Keisler, C. L. Reichardt, K. A. Aird, B. A. Benson, 
L. E. Bleem, J. E. Carlstrom, C. L. Chang, H. M. Cho, 

T. M. Crawford, A. T. Crites, T. de Haan, M. A. Dobbs, 

J. Dudley, E. M. George, N. W. Halverson, G. P. Holder, 
W. L. Holzapfel, S. Hoover, Z. Hou, J. D. Hrubes, M. Joy, 

L. Knox, A. T. Lee, E. M. Leitch, M. Lueker, D. Luong- 
Van, J. J. McMahon, J. Mehl, S. S. Meyer, M. Mil- 
lea, J. J. Mohr, T. E. Montroy, T. Natoli, S. Padin, 
T. Plagge, C. Pryke, J. E. Ruhl, K. K. Schaffer, L. Shaw, 
E. Shirokoff, H. G. Spieler, Z. Staniszewski, A. A. Stark, 

K. Story, A. van Engelen, K. Vanderlinde, J. D. Vieira, 

R. Williamson, and O. Zahn, The Astrophysical Journal 
743, 28 (2011). 

[24] S. Das, T. A. Marriage, P. A. R. Ade, P. Aguirre, 

M. Amiri, J. W. Appel, L. F. Barrientos, E. S. Battistelli, 
J. R. Bond, B. Brown, B. Burger, J. Chervenak, M. J. 
Devlin, S. R. Dicker, W. Bertrand Doriese, J. Dunk- 
ley, R. Drinner, T. Essinger-Hileman, R. P. Fisher, 
J. W. Fowler, A. Hajian, M. Halpern, M. Hasselfield, 

C. Hernandez-Monteagudo, G. C. Hilton, M. Hilton, 

A. D. Hincks, R. Hlozek, K. M. Huffenberger, D. H. 
Hughes, J. P. Hughes, L. Infante, K. D. Irwin, J. Bap¬ 
tiste Juin, M. Kaul, J. Klein, A. Kosowsky, J. M. Lau, 
M. Limon, Y.-T. Lin, R. H. Lupton, D. Marsden, K. Mar- 
tocci, P. Mauskopf, F. Menanteau, K. Moodley, H. Mose¬ 
ley, C. B. Netterfield, M. D. Niemack, M. R. Nolta, L. A. 
Page, L. Parker, B. Partridge, B. Reid, N. Sehgal, B. D. 
Sherwin, J. Sievers, D. N. Spergel, S. T. Staggs, D. S. 
Swetz, E. R. Switzer, R. Thornton, H. Trac, C. Tucker, 

R. Warne, E. Wollack, and Y. Zhao, Astrophys. J. 729, 
62 (2011). 

[25] K. T. Story, C. L. Reichardt, Z. Hou, R. Keisler, K. A. 
Aird, B. A. Benson, L. E. Bleem, J. E. Carlstrom, 

C. L. Chang, H.-M. Cho, T. M. Crawford, A. T. Crites, 
T. de Haan, M. A. Dobbs, J. Dudley, B. Follin, E. M. 
George, N. W. Halverson, G. P. Holder, W. L. Holzapfel, 

S. Hoover, J. D. Hrubes, M. Joy, L. Knox, A. T. Lee, 
E. M. Leitch, M. Lueker, D. Luong-Van, J. J. McMa¬ 
hon, J. Mehl, S. S. Meyer, M. Millea, J. J. Mohr, T. E. 
Montroy, S. Padin, T. Plagge, C. Pryke, J. E. Ruhl, 
J. T. Sayre, K. K. Schaffer, L. Shaw, E. Shirokoff, H. G. 
Spieler, Z. Staniszewski, A. A. Stark, A. van Enge¬ 
len, K. Vanderlinde, J. D. Vieira, R. Williamson, and 

O. Zahn, The Astrophysical Journal 779, 86 (2013). 

[26] S. Das, T. Louis, M. R. Nolta, G. E. Addison, E. S. 
Battistelli, J. R. Bond, E. Calabrese, D. Crichton, 
M. J. Devlin, S. Dicker, J. Dunkley, R. Diinner, J. W. 
Fowler, M. Gralla, A. Hajian, M. Halpern, M. Hasselfield, 
M. Hilton, A. D. Hincks, R. Hlozek, K. M. Huffenberger, 
J. P. Hughes, K. D. Irwin, A. Kosowsky, R. H. Lupton, 

T. A. Marriage, D. Marsden, F. Menanteau, K. Moodley, 
M. D. Niemack, L. A. Page, B. Partridge, E. D. Reese, 

B. L. Schmitt, N. Sehgal, B. D. Sherwin, J. L. Sievers, 

D. N. Spergel, S. T. Staggs, D. S. Swetz, E. R. Switzer, 
R. Thornton, H. Trac, and E. Wollack, Journal of Cos¬ 
mology and Astro-Particle Physics 04, 014 (2014). 

[27] P. Collaboration, P. A. R. Ade, N. Aghanim, 

C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio- 
Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, 


17 


R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, 

A. Benoit, A. Benoit-Levy, J.-P. Bernard, M. Bersanelli, 
P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi, J. R. 
Bond, J. Borrill, F. R. Bouchet, M. Bridges, M. Bucher, 

C. Burigana, R. C. Butler, E. Calabrese, B. Cappellini, 
J.-F. Cardoso, A. Catalano, A. Challinor, A. Chamballu, 

R. -R. Chary, X. Chen, H. C. Chiang, L.-Y. Chiang, P. R. 
Christensen, S. Church, D. L. Clements, S. Colombi, 

L. P. L. Colombo, F. Couc.hot, A. Coulais, B. P. Crill, 
A. Curto, F. Cuttaia, L. Danese, R. D. Davies, R. J. 
Davis, P. de Bernardis, A. de Rosa, G. de Zotti, J. De- 
labrouille, J.-M. Delouis, F.-X. Desert, C. Dickinson, 
J. M. Diego, K. Dolag, H. Dole, S. Donzelli, O. Dore, 

M. Douspis, J. Dunkley, X. Dupac, G. Efstathiou, F. Ei¬ 
sner, T. A. Enfilin, H. K. Eriksen, F. Finelli, O. Forni, 
M. Frailis, A. A. Fraisse, E. Franceschi, T. C. Gaier, 

S. Galeotta, S. Galli, K. Ganga, M. Giard, G. Gi- 
ardino, Y. Giraud-Heraud, E. Gjerlpw, J. Gonzalez- 
Nuevo, K. M. Gorski, S. Gratton, A. Gregorio, A. Grup- 
puso, J. E. Gudmundsson, J. Haissinski, J. Hamann, 
F. K. Hansen, D. Hanson, D. Harrison, S. Henrot- 
Versille, C. Hernandez-Monteagudo, D. Herranz, S. R. 
Hildebrandt, E. Hivon, M. Hobson, W. A. Holmes, 
A. Hornstrup, Z. Hou, W. Hovest, K. M. Huffenberger, 

A. H. Jaffe, T. R. Jaffe, J. Jewell, W. C. Jones, M. Ju- 
vela, E. Keihanen, R. Keskitalo, T. S. Kisner, R. Kneissl, 
J. Knoche, L. Knox, M. Kunz, H. Kurki-Suonio, G. La- 
gache, A. Lahteenmaki, J.-M. Lamarre, A. Lasenby, 
M. Lattanzi, R. J. Laureijs, C. R. Lawrence, S. Leach, 
J. P. Leahy, R. Leonardi, J. Leon-Tavares, J. Lesgour- 
gues, A. Lewis, M. Liguori, P. B. Lilje, M. Linden-Vprnle, 
M. Lopez-Caniego, P. M. Lubin, J. F. Macfas-Perez, 

B. Maffei, D. Maino, N. Mandolesi, M. Maris, D. J. 
Marshall, P. G. Martin, E. Martinez-Gonzalez, S. Masi, 
M. Massardi, S. Matarrese, F. Matthai, P. Mazzotta, 
P. R. Meinhold, A. Melchiorri, J.-B. Melin, L. Mendes, 

E. Menegoni, A. Mennella, M. Migliaccio, M. Millea, 

S. Mitra, M.-A. Miville-Deschenes, A. Moneti, L. Mon¬ 
tier, G. Morgante, D. Mortlock, A. Moss, D. Munshi, 
J. A. Murphy, P. Naselsky, F. Nati, P. Natoli, C. B. Net- 
terheld, H. U. Nprgaard-Nielsen, F. Noviello, D. Novikov, 

I. Novikov, I. J. O’Dwyer, S. Osborne, C. A. Oxbor- 
row, F. Paci, L. Pagano, F. Pajot, D. Paoletti, B. Par¬ 
tridge, F. Pasian, G. Patanchon, D. Pearson, T. J. Pear¬ 
son, H. V. Peiris, O. Perdereau, L. Perotto, F. Per- 
rotta, V. Pettorino, F. Piacentini, M. Piat, E. Pierpaoli, 

D. Pietrobon, S. Plaszczynski, P. Platania, E. Pointe- 
couteau, G. Polenta, N. Ponthieu, L. Popa, T. Pouta- 
nen, G. W. Pratt, G. Prezeau, S. Prunet, J.-L. Puget, 

J. P. Rachen, W. T. Reach, R. Rebolo, M. Reinecke, 
M. Remazeilles, C. Renault, S. Ricciardi, T. Riller, I. Ris- 
torcelli, G. Rocha, C. Rosset, G. Roudier, M. Rowan- 
Robinson, J. A. Rubino Martin, B. Rushohne, M. San- 
dri, D. Santos, M. Savelainen, G. Savini, D. Scott, M. D. 
Seiffert, E. P. S. Shellard, L. D. Spencer, J.-L. Starck, 
V. Stolyarov, R. Stompor, R. Sudiwala, R. Sunyaev, 

F. Sureau, D. Sutton, A.-S. Suur-Uski, J.-F. Sygnet, 
J. A. Tauber, D. Tavagnacco, L. Terenzi, L. Toffo- 
latti, M. Tomasi, M. Tristram, M. Tucci, J. Tuovi- 
nen, M. Tiirler, G. Umana, L. Valenziano, J. Valiviita, 
B. Van Tent, P. Vielva, F. Villa, N. Vittorio, L. A. Wade, 
B. D. Wandelt, I. K. Wehus, M. White, S. D. M. White, 
A. Wilkinson, D. Yvon, A. Zacchei, and A. Zonca, ArXiv 
e-prints 1303, 5076 (2013). 


[28] P. collaboration, P. A. R. Ade, N. Aghanim, C. Armitage- 
Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, 
J. Aumont, C. Baccigalupi, A. J. Banday, R. B. 

Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, 

A. Benoit, A. Benoit-Levy, J.-P. Bernard, M. Bersanelli, 
P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi, 

L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, 

F. Boulanger, M. Bridges, M. Bucher, C. Burigana, 

R. C. Butler, E. Calabrese, J.-F. Cardoso, A. Cata¬ 
lano, A. Challinor, A. Chamballu, L.-Y. Chiang, H. C. 
Chiang, P. R. Christensen, S. Church, D. L. Clements, 

S. Colombi, L. P. L. Colombo, C. Combet, F. Couchot, 

A. Coulais, B. P. Crill, A. Curto, F. Cuttaia, L. Danese, 

R. D. Davies, R. J. Davis, P. de Bernardis, A. de Rosa, 

G. de Zotti, J. Delabrouille, J.-M. Delouis, F.-X. Desert, 

C. Dickinson, J. M. Diego, H. Dole, S. Donzelli, O. Dore, 

M. Douspis, J. Dunkley, X. Dupac, G. Efstathiou, F. Ei¬ 
sner, T. A. Ensslin, H. K. Eriksen, F. Finelli, O. Forni, 
M. Frailis, A. A. Fraisse, E. Franceschi, T. C. Gaier, 

S. Galeotta, S. Galli, K. Ganga, M. Giard, G. Gi- 
ardino, Y. Giraud-Heraud, E. Gjerlow, J. Gonzalez- 
Nuevo, K. M. Gorski, S. Gratton, A. Gregorio, A. Grup- 
puso, J. E. Gudmundsson, F. K. Hansen, D. Hanson, 

D. Harrison, G. Helou, S. Henrot-Versille, C. Hernandez- 
Monteagudo, D. Herranz, S. R. Hildebrandt, E. Hivon, 
M. Hobson, W. A. Holmes, A. Hornstrup, W. Hov¬ 
est, K. M. Huffenberger, G. Hurier, T. R. Jaffe, A. H. 
Jaffe, J. Jewell, W. C. Jones, M. Juvela, E. Keiha¬ 
nen, R. Keskitalo, K. Kiiveri, T. S. Kisner, R. Kneissl, 
J. Knoche, L. Knox, M. Kunz, H. Kurki-Suonio, G. La- 
gache, A. Lahteenmaki, J.-M. Lamarre, A. Lasenby, 
M. Lattanzi, R. J. Laureijs, C. R. Lawrence, M. L. 
Jeune, S. Leach, J. P. Leahy, R. Leonardi, J. Leon- 
Tavares, J. Lesgourgues, M. Liguori, P. B. Lilje, V. Lind- 
holm, M. Linden-Vornle, M. Lopez-Caniego, P. M. Lubin, 
J. F. Macias-Perez, B. Maffei, D. Maino, N. Mandolesi, 

D. Marinucci, M. Maris, D. J. Marshall, P. G. Martin, 

E. Martinez-Gonzalez, S. Masi, S. Matarrese, F. Matthai, 
P. Mazzotta, P. R. Meinhold, A. Melchiorri, L. Mendes, 

E. Menegoni, A. Mennella, M. Migliaccio, M. Millea, 

S. Mitra, M.-A. Miville-Deschenes, D. Molinari, A. Mon¬ 
eti, L. Montier, G. Morgante, D. Mortlock, A. Moss, 
D. Munshi, P. Naselsky, F. Nati, P. Natoli, C. B. Netter- 
field, H. U. Norgaard-Nielsen, F. Noviello, D. Novikov, 

I. Novikov, I. J. O’Dwyer, F. Orieux, S. Osborne, C. A. 
Oxborrow, F. Paci, L. Pagano, F. Pajot, R. Paladini, 
D. Paoletti, B. Partridge, F. Pasian, G. Patanchon, 
P. Paykari, O. Perdereau, L. Perotto, F. Perrotta, F. Pia¬ 
centini, M. Piat, E. Pierpaoli, D. Pietrobon, S. Plaszczyn¬ 
ski, E. Pointecouteau, G. Polenta, N. Ponthieu, L. Popa, 

T. Poutanen, G. W. Pratt, G. Prezeau, S. Prunet, 

J. -L. Puget, J. P. Rachen, A. Rahlin, R. Rebolo, 
M. Reinecke, M. Remazeilles, C. Renault, S. Ricciardi, 

T. Riller, C. Ringeval, I. Ristorcelli, G. Rocha, C. Ros¬ 
set, G. Roudier, M. Rowan-Robinson, J. A. Rubino- 
Martin, B. Rusholme, M. Sandri, L. Sanselme, D. San¬ 
tos, G. Savini, D. Scott, M. D. Seiffert, E. P. S. Shellard, 
L. D. Spencer, J.-L. Starck, V. Stolyarov, R. Stompor, 
R. Sudiwala, F. Sureau, D. Sutton, A.-S. Suur-Uski, J.-F. 
Sygnet, J. A. Tauber, D. Tavagnacco, L. Terenzi, L. Tof- 
folatti, M. Tomasi, M. Tristram, M. Tucci, ,J. Tuovi- 
nen, M. Turler, L. Valenziano, J. Valiviita, B. Van Tent, 
J. Varis, P. Vielva, F. Villa, N. Vittorio, L. A. Wade, 

B. D. Wandelt, I. K. Wehus, M. White, S. D. M. White, 


18 


D. Yvon, A. Zacchei, and A. Zonca, arXiv: 1303.5075 
[astro-ph] (2013). 

[29] E. Aver, K. A. Olive, R. L. Porter, and E. D. Skillman, 
Journal of Cosmology and Astro-Particle Physics 11, 017 
(2013). 

[30] F. Iocco, G. Mangano, G. Miele, O. Pisanti, and P. D. 
Serpico, Physics Reports 472, 1 (2009), arXiv:0809.0631. 

[31] G. Steigman, Annual Review of Nuclear and Particle Sci¬ 
ence 57, 463 (2007), arXiv:0712.1100. 

[32] B. D. Fields, Annual Review of Nuclear and Particle Sci¬ 
ence 61, 47 (2011), arXiv: 1203.3551. 

[33] M. Pospelov and J. Pradler, Annual Review of Nuclear 
and Particle Science 60, 539 (2010). 

[34] A. Arbey, Computer Physics Communications 183, 1822 
(2012), arXiv: 1106.1363. 

[35] R. H. Cyburt, J. Ellis, B. D. Fields, F. Luo, K. A. Olive, 
and V. C. Spanos, Journal of Cosmology and Astro- 
Particle Physics 10, 021 (2009). 

[36] R. J. Cooke, M. Pettini, R. A. Jorgenson, M. T. Murphy, 
and C. C. Steidel, Astrophys. J. 781, 31 (2014). 

[37] G. Fiorentini, E. Lisi, S. Sarkar, and F. L. Villante, Phys¬ 
ical Review D 58 (1998), 10.1103/PhysRevD.58.063506, 
arXiv: astro-ph/9803177. 

[38] Y. I. Izotov and T. X. Thuan, The Astrophysical Journal 
Letters 710, L67 (2010). 

[39] F. Spite and M. Spite, Astron. & Astrophys. 115, 357 
(1982). 

[40] L. Sbordone, P. Bonifacio, E. Caffau, H.-G. Ludwig, 

N. T. Behara, J. I. Gonzalez Hernandez, M. Steffen, 

R. Cayrel, B. Freytag, C. van’t Veer, P. Molaro, B. Plez, 
T. Sivarani, M. Spite, F. Spite, T. C. Beers, N. Christlieb, 
P. FranCcois, and V. Hill, Astron. & Astrophys. 522, 
A26 (2010), arXiv:1003.4510 [astro-ph.GA], 

[41] J. Ellis, J. E. Kim, and D. V. Nanopoulos, Physics Let¬ 
ters B 145, 181 (1984). 

[42] M. Kawasaki and T. Moroi, Progress of Theoretical 
Physics 93, 879 (1995), hep-ph/9403364. 

[43] R. H. Cyburt, J. Ellis, B. D. Fields, and K. A. Olive, 
Phys. Rev. D 67. 103521 (2003), astro-ph/0211258. 

[44] M. Kleban and R. Rabadan, ArXiv High Energy Physics 
- Phenomenology e-prints (2005). 

[45] J. D. Bjorken, S. Ecklund, W. R. Nelson, A. Abashian, 
C. Church, B. Lu, L. W. Mo, T. A. Nunamaker, and 
P. Rassmann, Phys.Rev. D38, 3375 (1988). 

[46] G. G. Raffelt and D. S. P. Dearborn, Physical Review D 
37, 549 (1988). 

[47] G. G. Raffelt, Stars as Laboratories for Fundamen¬ 
tal Physics: The Astrophysics of Neutrinos, Axions, 
and Other Weakly Interacting Particles (University of 
Chicago Press, Chicago, 1996). 

[48] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979). 

[49] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, 
Nuclear Physics B 166, 493 (1980). 

[50] A. R. Zhitnitskij, (1980). 

[51] M. Dine, W. Fischler, and M. Srednicki, Physics Letters 
B 104, 199 (1981). 

[52] J. Beringer, J.-F. Arguin, R. M. Barnett, K. Copic, 

O. Dahl, D. E. Groom, C.-J. Lin, J. Lys, H. Mu- 
rayama, C. G. Wohl, W.-M. Yao, P. A. Zyla, C. Am- 
sler, M. Antonelli, D. M. Asner, H. Baer, H. R. Band, 
T. Basaglia, C. W. Bauer, J. J. Beatty, V. I. Be¬ 
lousov, E. Bergren, G. Bernardi, W. Bertl, S. Bethke, 
H. Bichsel, O. Biebel, E. Blucher, S. Blusk, G. Brooij- 
mans, O. Buchmueller, R. N. Cahn, M. Carena, A. Cec- 


cucci, D. Chakraborty, M.-C. Chen, R. S. Chivukula, 
G. Cowan, G. D’Ambrosio, T. Damour, D. de Florian, 

A. de Gouvea, T. DeGrand, P. de Jong, G. Dissertori, 

B. Dobrescu, M. Doser, M. Drees, D. A. Edwards, S. Ei- 
delman, J. Erler, V. V. Ezhela, W. Fetscher, B. D. 
Fields, B. Foster, T. K. Gaisser, L. Garren, H.-J. Gerber, 
G. Gerbier, T. Gherghetta, S. Golwala, M. Goodman, 

C. Grab, A. V. Gritsan, J.-F. Grivaz, M. Griinewald, 
A. Gurtu, T. Gutsche, H. E. Haber, K. Hagiwara, 

C. Hagmann, C. Hanhart, S. Hashimoto, K. G. Hayes, 
M. Heffner, B. Heltsley, J. J. Hernandez-Rey, K. Hikasa, 
A. Hocker, J. Holder, A. Holtkamp, J. Huston, J. D. 
Jackson, K. F. Johnson, T. Junk, D. Karlen, D. Kirkby, 

S. R. Klein, E. Klempt, R. V. Kowalewski, F. Krauss, 
M. Kreps, B. Krusche, Y. V. Kuyanov, Y. Kwon, O. La- 
hav, J. Laiho, P. Langacker, A. Liddle, Z. Ligeti, T. M. 
Liss, L. Littenberg, K. S. Lugovsky, S. B. Lugovsky, 

T. Mannel, A. V. Manohar, W. J. Marciano, A. D. Mar¬ 
tin, A. Masoni, J. Matthews, D. Milstead, R. Miquel, 
K. Monig, F. Moortgat, K. Nakamura, M. Narain, P. Na¬ 
son, S. Navas, M. Neubert, P. Nevski, Y. Nir, K. A. 
Olive, L. Pape, J. Parsons, C. Patrignani, J. A. Pea¬ 
cock, S. T. Petcov, A. Piepke, A. Pomarol, G. Punzi, 
A. Quadt, S. Raby, G. Raffelt, B. N. Ratcliff, P. Richard¬ 
son, S. Roesler, S. Rolli, A. Romaniouk, L. J. Rosen¬ 
berg, J. L. Rosner, C. T. Sachrajda, Y. Sakai, G. P. 
Salam, S. Sarkar, F. Sauli, O. Schneider, K. Scholberg, 

D. Scott, W. G. Seligman, M. H. Shaevitz, S. R. Sharpe, 
M. Silari, T. Sjostrand, P. Skands, J. G. Smith, G. F. 
Smoot, S. Spanier, H. Spieler, A. Stahl, T. Stanev, 

S. L. Stone, T. Sumiyoshi, M. J. Syphers, F. Takahashi, 

M. Tanabashi, J. Terning, M. Titov, N. P. Tkachenko, 

N. A. Tornqvist, D. Tovey, G. Valencia, K. van Bib¬ 
ber, G. Venanzoni, M. G. Vincter, P. Vogel, A. Vogt, 
W. Walkowiak, C. W. Walter, D. R. Ward, T. Watari, 
G. Weiglein, E. J. Weinberg, L. R. Wiencke, L. Wolfen- 
stein, J. Womersley, C. L. Woody, R. L. Workman, A. Ya¬ 
mamoto, G. P. Zeller, O. V. Zenin, J. Zhang, R.-Y. Zhu, 

G. Harper, V. S. Lugovsky, and P. Schaffner, Phys. Rev. 
D 86, 010001 (2012). 

[53] D. B. Kaplan, Nuclear Physics B 260, 215 (1985). 

[54] S. L. Cheng, C. Q. Geng,' and W.-T. Ni, Phys. Rev. D 
52, 3132 (1995). 

[55] https : //github. com/marius311/cosmoslik. 

[56] C. M. Ho and R. J. Scherrer, Phys. Rev. D 87, 065016 
(2013). 

[57] K. N. Abazajian, K. Arnold, J. Austermann, B. A. 
Benson, C. Bischoff, J. Bock, J. R. Bond, J. Bor- 
rill, E. Calabrese, J. E. Carlstrom, C. S. Carvalho, 
C. L. Chang, H. C. Chiang, S. Church, A. Cooray, 

T. M. Crawford, K. S. Dawson, S. Das, M. J. Devlin, 
M. Dobbs, S. Dodelson, O. Dore, J. Dunkley, J. Errard, 
A. Fraisse, J. Gallicchio, N. W. Halverson, S. Hanany, 
S. R. Hildebrandt, A. Hincks, R. Hlozek, G. Holder, 
W. L. Holzapfel, K. Honscheid, W. Hu, J. Hubmayr, 
K. Irwin, W. C. Jones, M. Kamionkowski, B. Keating, 

R. Keisler, L. Knox, E. Komatsu, J. Kovac, C.-L. Kuo, 
C. Lawrence, A. T. Lee, E. Leitch, E. Linder, P. Lubin, 
J. McMahon, A. Miller, L. Newburgh, M. D. Niemack, 

H. Nguyen, H. T. Nguyen, L. Page, C. Pryke, C. L. 
Reichardt, J. E. Ruhl, N. Sehgal, U. Seljak, J. Sievers, 

E. Silverstein, A. Slosar, K. M. Smith, D. Spergel, S. T. 
Staggs, A. Stark, R. Stompor, A. G. Vieregg, G. Wang, 

S. Watson, E. J. Wollack, W. L. K. Wu, K. W. Yoon, 


19 


and O. Zahn, ArXiv e-prints 1309, 5383 (2013). 

[58] Their forecast assumes standard BBN, however we be¬ 
lieve it unlikely that freeing Y p could degrade it so much 
as to allow 2.44. 

[59] T. B. Collaboration, Physical Review Letters 113 (2014), 
10.1103/PhysRevLett.113.201801, arXiv: 1406.2980. 

[60] T. Abe, I. Adachi, K. Adamczyk, S. Ahn, H. Aihara, 

K. Akai, M. Aloi, L. Andricek, K. Aoki, Y. Arai, 

A. Arefiev, K. Arinstein, Y. Arita, D. M. Asner, 
V. Aulchenko, T. Aushev, T. Aziz, A. M. Bakich, V. Bal- 
agura, Y. Ban, E. Barberio, T. Barvich, K. Belous, 
T. Bergauer, V. Bhardwaj, B. Bhuyan, S. Blyth, A. Bon¬ 
dar, G. Bonvicini, A. Bozek, M. Bracko, J. Brodz- 
icka, O. Brovchenko, T. E. Browder, G. Cao, M.-C. 
Chang, P. Chang, Y. Chao, V. Chekelian, A. Chen, K.-F. 
Chen, P. Chen, B. G. Cheon, C.-C. Chiang, R. Chistov, 

K. Cho, S.-K. Choi, K. Chung, A. Comerma, M. Cooney, 
D. E. Cowley, T. Critchlow, J. Dalseno, M. Danilov, 

A. Dieguez, A. Dierlamm, M. Dillon, J. Dingfelder, 

R. Dolenec, Z. Dolezal, Z. Drasal, A. Drutskoy, W. Dun- 
gel, D. Dutta, S. Eidelman, A. Enomoto, D. Epifanov, 

S. Esen, J. E. Fast, M. Feindt, M. F. Garcia, T. Fi- 
held, P. Fischer, J. Flanagan, S. Fourletov, J. Fourletova, 

L. Freixas, A. Frey, M. Friedl, R. Fruehwirth, H. Fujii, 

M. Fujikawa, Y. Fukuma, Y. Funakoshi, K. Furukawa, 
J. Fuster, N. Gabyshev, A. G. d. V. Cueto, A. Gar¬ 
mash, L. Garrido, C. Geisler, I. Gfall, Y. M. Goh, 

B. Golob, I. Gorton, R. Grzymkowski, H. Guo, H. Ha, 

J. Haba, K. Hara, T. Hara, T. Haruyama, K. Hayasaka, 

K. Hayashi, H. Hayashii, M. Heck, S. Heindl, C. Heller, 

T. Hemperek, T. Higuchi, Y. Horii, W.-S. Hou, Y. B. Hsi- 
ung, C.-H. Huang, S. Hwang, H. J. Hyun, Y. Igarashi, 

C. Iglesias, Y. Iida, T. Iijima, M. Imamura, K. Inami, 
C. Irmler, M. Ishizuka, K. Itagaki, R. Itoh, M. Iwabuchi, 

G. Iwai, M. Iwai, M. Iwasaki, M. Iwasaki, Y. Iwasaki, 
T. Iwashita, S. Iwata, H. Jang, X. Ji, T. Jinno, M. Jones, 
T. Julius, T. Kageyama, D. H. Kah, H. Kakuno, 
T. Kamitani, K. Kanazawa, P. Kapusta, S. U. Kataoka, 

N. Katayama, M. Kawai, Y. Kawai, T. Kawasaki, 

J. Kennedy, H. Kichimi, M. Kikuchi, C. Kiesling, B. K. 
Kim, G. N. Kim, H. J. Kim, H. O. Kim, J.-B. Kim, J. H. 
Kim, M. J. Kim, S. K. Kim, K. T. Kim, T. Y. Kim, 

K. Kinoshita, K. Kishi, B. Kisielewski, K. K. van Dam, 
J. Knopf, B. R. Ko, M. Koch, P. Kodys, C. ICoffmane, 
Y. Koga, T. Kohriki, S. Koike, H. Koiso, Y. Kondo, 
S. Korpar, R. T. Kouzes, C. Kreidl, M. Kreps, P. Krizan, 
P. Krokovny, H. Krueger, A. Kruth, W. Kuhn, T. Kuhr, 

R. Kumar, T. Kumita, S. Kupper, A. Kuzmin, P. Kvas- 
nicka, Y.-J. Kwon, C. Lacasta, J. S. Lange, I.-S. Lee, 
M. J. Lee, M. W. Lee, S.-H. Lee, M. Lemarenko, J. Li, 


W. D. Li, Y. Li, J. Libby, A. Limosani, C. Liu, H. Liu, 

Y. Liu, Z. Liu, D. Liventsev, A. L. Virto, Y. Makida, 

Z. P. Mao, C. Marinas, M. Masuzawa, D. Matvienko, 
W. Mitaroff, K. Miyabayashi, H. Miyata, Y. Miyazaki, 
T. Miyoshi, R. Mizuk, G. B. Mohanty, D. Mohap- 
atra, A. Moll, T. Mori, A. Morita, Y. Morita, H.- 

G. Moser, D. M. Martin, T. Mueller, D. Muenchow, 

J. Murakami, S. S. Myung, T. Nagamine, I. Nakamura, 
T. T. Nakamura, E. Nakano, H. Nakano, M. Nakao, 

H. Nakazawa, S.-H. Nam, Z. Natkaniec, E. Nedelkovska, 

K. Negishi, S. Neubauer, C. Ng, J. Ninkovic, S. Nishida, 
K. Nishimura, E. Novikov, T. Nozaki, S. Ogawa, 
K. Ohmi, Y. Ohnishi, T. Ohshima, N. Ohuchi, K. Oide, 

S. L. Olsen, M. Ono, Y. Ono, Y. Onuki, W. Ostrowicz, 
H. Ozaki, P. Pakhlov, G. Pakhlova, H. Palka, H. Park, 
H. K. Park, L. S. Peak, T. Peng, I. Peric, M. Pernicka, 

R. Pestotnik, M. Petrie, L. E. Piilonen, A. Poluektov, 
M. Prim, K. Prothmann, K. Regimbal, B. Reisert, R. H. 
Richter, J. Riera-Babures, A. Ritter, A. Ritter, M. Ritter, 

M. Roehrken, J. Rorie, M. Rosen, M. Rozanska, L. Ruck- 
man, S. Rummel, V. Rusinov, R. M. Russell, S. Ryu, 
H. Sahoo, K. Sakai, Y. Sakai, L. Santelj, T. Sasaki, 

N. Sato, Y. Sato, J. Scheirich, J. Schieck, C. Schwanda, 
A. J. Schwartz, B. Schwenker, A. Seljak, K. Senyo, O.- 

S. Seon, M. E. Sevior, M. Shapkin, V. Shebalin, C. P. 
Shen, H. Shibuya, S. Shiizuka, J.-G. Shiu, B. Shwartz, 

F. Simon, H. J. Simonis, J. B. Singh, R. Sinha, M. Sitarz, 
P. Smerkol, A. Sokolov, E. Solovieva, S. Stanic, M. Staric, 

J. Stypula, Y. Suetsugu, S. Sugihara, T. Sugimura, 

K. Sumisawa, T. Sumiyoshi, K. Suzuki, S. Y. Suzuki, 

H. Takagaki, F. Takasaki, H. Takeichi, Y. Takubo, 
M. Tanaka, S. Tanaka, N. Taniguchi, E. Tarkovsky, 

G. Tatishvili, M. Tawada, G. N. Taylor, Y. Teramoto, 

I. Tikhomirov, K. Trabelsi, T. Tsuboyama, K. Tsunada, 
Y.-C. Tu, T. Uchida, S. Uehara, K. Ueno, T. Uglov, 

Y. Unno, S. Uno, P. Urquijo, Y. Ushiroda, Y. Usov, 

S. Vahsen, M. Valentan, P. Vanhoefer, G. Varner, K. E. 
Varvell, P. Vazquez, I. Vila, E. Vilella, A. Vinokurova, 

J. Visniakov, M. Vos, C. H. Wang, J. Wang, M.-Z. 
Wang, P. Wang, A. Wassatch, M. Watanabe, Y. Watase, 

T. Weiler, N. Wermes, R. E. Wescott, E. White, J. Wicht, 

L. Widhalm, K. M. Williams, E. Won, H. Xu, B. D. Yab- 
sley, H. Yamamoto, H. Yamaoka, Y. Yamaoka, M. Ya- 
mauchi, Y. Yin, H. Yoon, J. Yu, C. Z. Yuan, Y. Yusa, 

D. Zander, M. Zdybal, Z. P. Zhang, J. Zhao, L. Zhao, 

Z. Zhao, V. Zhilich, P. Zhou, V. Zhulanov, T. Zivko, 
A. Zupanc, and O. Zyukova, arXiv:1011.0352 [hep-ex, 
physics:physics] (2010), arXiv: 1011.0352. 



